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This review presents the mathematical aspect of the work done on level cross- 
ings and upcrossings of a level, and insists on the different approaches used to 
tackle a given problem. We speak only briefly about examples or applications 
which of course were the reason for such theoretical work and which have been 
greatly developing in the extreme value field over the past years. This review 
of mathematical results and methods can also be very useful in mathematical 
statistics applications, reliability theory and other areas of applications of 
stochastic processes. 

When possible, we choose to privilege conditions on the behavior of spectral 
or covariance functions rather than general conditions of mixing, since they are 
easier to work with. 

We mainly consider continuous parameter processes, since they appear in most 
of the mathematical models that describe physical phenomena. Nevertheless one 
approach to work on level crossings or on upcrossings of a level that are analo- 
gous to the exceedances used in the discrete case, can be through discretization, 
using discrete parameter results. 

Let X — (Xt,t € M d ) be a real stochastic process. Our main concern is the 
measure of the random set of level x defined by '■= {t : X t — x}, by taking 
the same notation as Wschebor (see |163j ). The unidimensional case will be ex- 
posed in detail with results and methods (especially the ones which can/could 
be easily adapted to a higher dimension), whereas only a partial view of the 
multidimensional case will be given. 

As a preliminary, let us mention two important results (see [31] for the first 
one and for instance |90j for the second one) which reinforce our interest in the 
study of C*. 

Theorem 1.1 (Bulinskaya, 1961) 

Let X = (X s ,0 < s < t) be an a.s. continuously differentiable stochastic process 
with one- dimensional density p s (u) bounded in u for all s. 
Then for any level x the probability of existence of a point s such that the events 
(X(s) = x) and (X(s) = 0) occur simultaneously, is equal to 0. In particular, 
the probability that X t becomes tangent to the level x, is equal to 0, in other 
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words the probability of contingence of the level x by the process X is equal to 
0. 

In the case of a Gaussian process, as a consequence of the separability property 
and the Tsyrelson theorem, we have 

Theorem 1.2 For any Gaussian process X on a arbitrary parameter set T , for 
any x such that 

P[supX(s) < x]>0, 

sGT 

there is no contingence, with probability one. 

Let us give the notations and specific symbols used through this survey. 

L p (fl) := {X : E[\X\ P ] < oo}; ||.|| p : L^-norm. 

1a'- indicator function of a set A. 

4> Si t(x,x,y,y): density function of (X s , X s , X t , X t ); 

cf)t{x,x,y,y) := (f)o, t (x,x,y,y). 

r.v.: random variable 

CLT: central limit theorem 

s.t.: such that 

— > stands for convergence in distribution. 
2. Crossings of Gaussian processes 

Studies on level-crossings by stationary Gaussian processes began about sixty 
years ago. Different approaches have been proposed. Here is a survey of the 
literature on the number of crossings of a given level or of a differentiable curve 
in a fixed time interval by a continuous spectrum Gaussian process. 
Besides the well-known books about Extremes, let us quote a short survey by 
Slud in 1994 (see |148j or |149j ). a more general survey about extremes including 
a short section about level crossings by Leadbetter et al. in 1988 (see [55]) and 
by Rootzen in 1995 (see |132] L and another one by Piterbarg in his 1996 book 
(see [1 22] . in particular for some methods described in more detail than here). 
Our purpose here is to focus only on crossing counts in order to be more ex- 
plicit about the subject, not only recalling the main results, but also giving the 
main ideas about the different methods used to establish them. To make the 
methods easier to understand, we consider mainly the crossing counts of a given 
level. Note that the problem of curve (ip) crossings by a stationary process X 
may also be regarded as a zero-crossing problem for the non stationary process 
X* := X — tjj (but stationary in the sense of the covariance), as pointed out by 
Cramer and Leadbetter (see |38j). 

Let X = (X s ,s > 0) be a real stationary Gaussian process with variance one 

and a.s. continuous sample functions. 

Let tp be a continuously differentiable function. 

We denote Nj(x) or Ni(ip)) the number of crossings of a given level x or of a 



M.F. Kratz/Level crossings and other level Junctionals 



233 



curve respectively, X on the interval /, and Nf (x) the number of up- 
crossings of x by X (recall that X s is said to have an upcrossing of x at so > if 
for some e > 0, X s < x in (sq — e, sq) and X s > x in (sq, sq + e)). Hereafter this 
will be simply denoted N t (x) or N t (tp) or N^~(x) respectively, when I = [0,t], 
t > 0. 

£. J . Moments and factorial moments 

Distributional results about level or curve crossings by a Gaussian process are 
often obtained in terms of factorial moments for the number of crossings, that 
is why many authors have worked not only on moments but also on factorial 
moments to find out expressions or conditions (in terms of the covariance func- 
tion of the process) for their finiteness, since a certain number of applications 
require to know if they are finite. 

Note also that the conditions governing finite crossing moments are local ones, 
since Holder inequality implies that (E[N 2t (x)]) k < 2 k (E[N t (x)]) k , and there- 
fore when (E[N t (x)]) k is finite for some t, it is finite for all t. 

2.1.1. Introduction 

• Kac's method and formula (1943). 

We can start with Kac (see [73]), who studied the number N of zeros of a Gaus- 
sian random polynomial on some bounded interval of IR. To compute E[N] , he 
proposed a method which used, so to speak, in a formal way, the approximation 

of the "Dirac function" Sq(x) by the function — lr_ ££ )(x). He gave the first 

heuristic expression of N t (0) as a function of the process X: 

Theorem 2.1 (Kac formula, 1943) 
1 /°° /' 

N t (0) — — / / cos((X s )\X s \dsd£, with probability one. 

From now on, we suppose w.l.o.g. that X is centered, with correlation 
function r{t) = E[X X t ], given also in respect with the spectral function F as 

poo 

r{t) = I cos(\t)dF(X). We denote A2 the second moment (when it exists) of 
Jo 

/•OO 

the spectral function, i.e. A2 = / X 2 dF(\). 

Jo 

• Rice (1945). 

One of the best known first results is the one of Rice (see )129j ) who proved by 
intuitive methods related to those used later by Cramer and Leadbetter (as the 
discretization method described below), that 
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Theorem 2.2 (Rice formula, 1945) 



E[N t (x)] = te~ x 12 y/-r"(0) /tt. (2.1) 

It means that the mean number of crossings is the most important for zero 
crossings, and decreases exponentially with the level. 

• ltd (1964), Ylvisaker (1965). 

ltd and Ylvisaker have proved (|2.ip and therefore that the mean number of 
crossings is finite exactly when the second spectral moment is finite: 

Theorem 2.3 (ltd, 1964; Ylvisaker, 1965) 

E[N t (x)} < oo A 2 < oo <=> -r"(0) < oo. 

• The discretization method. 

From the intuitive method developed by Rice at the beginning of the 40s (see 
[129]), Ivanov in 1960 (see [70]), Bulinskaya in 1961 (see (3]]), Ito in 1964 (see 
[71] ) and Ylvisaker in 1965 (see [165] ) proposed rigorous proofs for zero count- 
ings. There followed the general formulation due to Leadbetter in 1966 (see [38] 
p. 195 or [85] p. 148), known as the method of discretization, which also applies 
to non-normal processes. 

This method is based on the approximation of the continuous process X — 
(X t , t > 0) by the sequence (X(jq),j = 1, 2, • • • ) where q satisfies some condi- 
tions related to the level x of upcrossings: 

q = q(x) — > and xq 1 clS X — X (t) — > oo (or equivalently as oo). 

(2.2) 

Consider any sequence (q n ) s.t. q n J. and the number N n j of points (jq n , j — 

U ~ l)3n S / and jq n G /, 
X((j'-l)g») <x<X(jq n ). 
It can be proved (see [38] p. 195, and [85], lemma 7.7.2) that 



1, 2, • • • ), in a fixed and bounded interval I, s.t. 



N n ,i — > Nj~(x) a.s. as n — > oo, (2.3) 

hence that Nf*~(x) is a (possibly infinite- valued) r.v., and that (see [85] ) 
JE[iV njI ] -» JE[iV+(x)]. 

n — ^oo 

Following the notation of Leadbetter et al. (see [85]), let 

J q (x) = -lP[X <x<X q }, q>0. (2.4) 

By choosing / = (0, t], we have v n := t/q n points jq n G I and by stationarity, it 
yields JE[N n ] — [y n — l)lP[Xo < X < X Qn ] ~ tJ qn (x); since the sequence (q n ) 
is arbitrary, it implies 

E[Nf(x)] =tlim J q (x) (2.5) 

£(J.O 
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and hence that 

E[N+{x)]=tE[N+{x)]. (2.6) 
Under mild conditions, lim J„(x), and therefore E[N^(x)], can be expressed 

in a simple integral form (by writing J q as P[Xq > x , X q > q^ 1 (x — Xq)], 
with X q := q~ 1 (X q — Xq), and by doing some change of variable), giving in the 
normal case (see [55], lemma 7.3.1): 

Lemma 2.1 (Leadbetter et al.) 

Let X = (X t ,t > 0) be a standardized stationary normal process with its second 
spectral moment A2 finite and let (q = q{x),x) satisfy (KM) - Then, as q — ► 0, 

J q {x) ~ / zp(x,z)dz = -|— e~ x2/2 , 
Jo 27T 

p(.,.) denoting the joint density of(X,X) and J q satisfying \2.^) . 

Therefore this last result, together with (|2.5|) , provide the Rice formula for the 
number of upcrossings of any fixed level x per unit time by a standardized 
stationary normal process, namely 

E[N+(x)} = e^ 2 / 2 A* /2 /(2^). 



• Cramer and Leadbetter (1965,1967), Ylvisaker (1966). 

In the 60s, following the work of Cramer, generalizations to curve crossings and 
higher order moments for N t (.) were considered in a series of papers by Cramer 
and Leadbetter (see [35]) and Ylvisaker (see |166j ). 
As regards curve crossings, the generalized Rice formula was obtained: 

Lemma 2.2 ( Generalized Rice formula; Ylvisaker, 1966; Cramer and Leadbet- 
ter, 1967). 

Let X = (X t , t > 0) be a mean zero, variance one, stationary Gaussian process, 
with twice differ entiable covariance function r. Let ip be a continuous differen- 
tiable real function. Then 

E[N t {i,)]=^f^)f y^{y)) 

J o 

where ip and $ are the standard normal density and distribution function 
respectively. 

Again the authors used the discretization method and approximated the contin- 
uous time number N t (ip) of ^-crossings by the discrete-time numbers (N^(t, q n )) 
of crossings of continuous polygonal curves agreeing with ip at points jq n 
(j = 0, 1, • • • , 2"), with time steps q n = t/2 n : 

2™-l 

Nj,(t,q n ):= 1 {(x ]qn ^{j qn ))(x (j+1)qn -i,([ 3 +i) qn ))<Q) ( 2 - 8 ) 




dy. 
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to obtain that N^(t, q n ) f N t (ip) with probability one, as n — > oo. □ 
2.1.2. Moments and factorial moments of order 2 

• Cramer and Leadbetter proposed a sufficient condition (known nowdays as the 
Geman condition) on the correlation function of X (stationary case) in order to 
have the random variable N t (x) belonging to L 2 (f2), namely 

Theorem 2.4 (Cramer et al, 1967) 

If 36 >0, L(t) := LMJZLM G £*([(), «5], die) (2.9) 
then E[N?(0)} < oo. 

• Explicit (factorial) moment formulas for the number of crossings have been 
obtained by Cramer and Leadbetter (see [35J ), Belyaev (see jT^]) and Ylvisaker 
(see [165] ) . based on careful computations involving joint densities of values and 
derivatives of the underlying Gaussian process. 

In particular, the second moment of N t (x) is given by 

E[N?(x)} = E[N t (x)}+2 f (t-u) f \x\\y\cj) u (x,x,x,y)dxdydu. 

Jo Jwe 

Concerning the second factorial moment, Cramer & Leadbetter provided an 
explicit formula for the number of zeros of the process X (see [35], p. 209), from 
which the following formula for the second factorial moment of the number of 
crossings of a continuous differentiable real function tp by X can be deduced: 

E[N t (i>)(N t (iP)-l)}=f f f Jii-^ tl ||i 3 -^t a |^i,t a (^ti > ii J V^.*a)d*id*a<fti«ft2 (2.10) 

Jo Jo Jwe 

where the density <j>t 1 ,t 2 is supposed non-singular for all t\ ^ t2- The formula 
holds whether the second factorial moment is finite or not. 



• Ershov (1967). 

This author proved (see [51]) that whenever a Gaussian stationary process X = 
(X t ,t S M) with covariance function r is with mixing (i.e. that \r(\i — j\)\ < 
f(\i — j\) with lim f(k) = 0), then the number of its x-upcrossings on all E 

k — >oo 

cannot hold a finite second moment: 

Theorem 2.5 (Ershov, 1967) 

If r(t) —>0 as t — > oo, 

then var (iV t + (x)) — > oo, as t — > oo. 



• Geman (1972). 

Geman proved in 1972 (see [S3]) that the condition (|2.9p was not only sufficient 
but also necessary, by showing that if r ^ diverges on (0, 5) then so 
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does the integral 

/OO poo 
/ \xy\(/>s(0, x, 0, y)dxdy appearing in the computation of the second 
- oo J — OO 

moment. Thus we have 

Theorem 2.6 (Geman condition, 1972) 

Geman condition (TQ|) ^ E[N?(0)] < oo. 

Of course, this result holds for any stationary Gaussian process when consider- 
ing the number of crossings of the mean of the process. 

• Wschebor (1985). 

This author provided in |163j (under the Geman condition) an explicit expres- 
sion in the case of two different levels x and y, with i/y, namely: 

E[N t (x)N t (y)} = (t-u) \x\\y\(tp u (x,x,y,y) + tp u (y,x,x,y))dxdydu. 
Jo JE 2 

Note that this expression differs from the one of Cramer et al. when x — > y, 
which means that the function lE[N t (x)N t (y)) is not continuous on the diagonal. 

• Piterbarg (1982), Kratz and Rootzen (1997). 

Let us assume that the correlation function r of X satisfies 

s 2 

r(s) = 1 - — + o(s 2 ) as s -> 0, (2.11) 

E{(X , (s)-X , (Q)) 2 } = 2(r , '(s)-r' , (0))<c\sp, s > 0, < 7 < 2, 

(2.12) 

that r(s) and its first derivative decay polynomially, 

\r(s)\ <Cs- a , \r(s)\+r'( S ) 2 <C QS - a , s>0, (2.13) 

for some a > 2 and constants c, C, Co, 
and that the range of t, x is such that 

£ < E[N+(t)] < K Q (2.14) 

for some fixed Kq, e > 0; by using Piterbarg's notations, let 



r'(s) i_ r(s )2_ r / (s)S 
rn = m{s) = f-^-, a = a\s) — 



l + r(s)' " " y "' y l-r(s) 2 
and 

./ \ m2 ( v ) * 0-2 ( s ) 

m (s) — sup =, a = sup- 



v>s>o 2ny/l — r 2 (v) s >q 2ny/l — r 2 (s) 

Piterbarg (see |121j ) proposed some bounds of the second factorial moment of 
the number of upcrossings; Kratz and Rootzen (see [M]) gave a small variation 
of his result, but with more precise bounds, namely 
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Lemma 2.3 (Kratz and Rootzen, 1997) 

Suppose that the previous hypotheses are satisfied (with 7 = 2). 
Then, for t, x > 1, 

E[N+(t)(N+(t) - 1)] < K 3 tx 2+2/a e-^( 1+ini ^ 0p{s) ) + K 4 t 2 x 2 e~ x2 , (2.15) 

with p{.) defined in [2~J2^ below, K 3 = (3.1a* + 2.3m*(0))(3C ) 1/c \ K 4 = 
2.6a* + 2.3to*(0). 

• Kratz and Leon (2006). 

Recently, Kratz and Leon (see [83]) have succeded in generalizing theorem l2.6l to 
any given level x (not only the mean of the process) and to some differentiable 
curve ip. 

Note that the problem of finding a simple necessary and sufficient condition for 
the number of crossings of any level has already been broached in some very 
interesting papers that proposed sufficient conditions as the ones of Cuzick (see 
[40] . [43], [11]). But getting necessary conditions remained an open problem 
for many years. At present, the solution of this problem is enunciated in the 
following theorem. 

Theorem 2.7 (Kratz and Leon, 2006) 

1 ) For any given level x, we have 

r"(t) - r"(0) 

E[N 2 (x)] < 00 & 36 > 0, L{t) = K ' ^ G L 1 ([0,S],dx) (Geman condition). 

2) Suppose that the continuous differentiable real function tp satisfies for some 
5>0, 

ds < 00, where j(t) is the modulus of continuity of ip. 
Then 

E[N?(if>)] < 00 & L(t) G L 1 ([0,S],dx). 

This smooth condition on ip is satisfied by a large class of functions which in- 
cludes in particular functions whose derivatives are Holder. 

The method used to prove that the Geman condition keeps being the sufficient 
and necessary condition to have a finite second moment is quite simple. 
It relies mainly on the study of some functions of r and its derivatives at the 
neighborhood of 0, and the chaos expansion of the second moment, a notion 
which will be explained later. 

2.1.3. Factorial moments and moments of higher order 

• Concerning moments of order higher than 2, Cramer and Leadbetter (see [38] ) 
got in the stationary case, under very mild conditions, results that Belyaev (see 
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|19j ) derived in the non-stationary case under slightly more restrictive condi- 
tions, weakened by Ylvisaker (see |166| ) in the stationary case but which may 
also be adapted to cover nonstationary cases. Let us give for instance the fcth 
factorial moment of the number of zero crossings by a stationary Gaussian pro- 
cess: 



M t fe (0):= E[N t (0)(N t (0) -!)■■• (N t (0) -k + 1)] 



Y[\£(tj)\ I *(*,■) = o, o<j<k 







[dh-- 


■ [dt k E 


la 


Jo 



p(0,...,0) (2.16) 



where p(x±, ■ • • ,Xk) is the joint density of the r.v. (X tl , • • • ,X tk ). 
• Belyaev (1967). 

Concurrently to this study, Belyaev proposed in |19j a sufficient condition for 
the finiteness of the fcth factorial moment Mj^(0) for the number N t (0) of zero 
crossings on the interval [0, t] in terms of the covariance matrix of (X tl , • • • , X tk ) 
and of 



var[X u | Xu =0, 1 <j < k 



(2.17) 



Theorem 2.8 (Belyaev, 1967) 



rt /tt* ^V 2 



If I dtf- I dt k I n ;-^ ) <oo (2.18) 



detT,k 

then M t fc (0) := ^iV t (0)(jV t (0) -!)•■• (N t (0) - k + 1)] < oo. 



• Cuzick (1975-1978). 

Cuzick proved in 1975 (see [30]) that Belyaev's condition (|2.18p for the finiteness 
of the fcth factorial moments for the number of zero crossings, was not only 
sufficient but also necessary: 

Theorem 2.9 (Cuzick condition, 1975) 

Condition (2lM [ dA t ■ ■ ■ [ dA k ^ ( Il^lf±] < oo (2.19) 

Jo Jo \ detY, k J 

for some e > 0, where Oi satisfies <\2. 171 
«*■ M t k (0) = E[N t (0){N t {0) - 1) ■ • • (JV t (0) - k + 1)] < oo. 

The proof that condition (|2.18[) is equivalent to condition (|2 . 1 9[) is immediate 
with the change of variables Aj := U + i — ti, 1 < i < k — 1, after noticing the 
symmetry of the integrand in (|2.18l) . 
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Now the necessity of (|2.19[) comes from the fact that the lemma given below im- 

r* t* (u k a 2 \ 1/2 

plies that ([STTfi]) > C dtf- dt k \ 1 , 1 1 , with C some constant. 
Jo Jo \ detY, k J 

Then Cuzick tried to derive from his result (|2.19p simpler sufficient conditions 
to have the finiteness of the /cth (factorial) moments for the number of crossings. 
In particular, in his 75s paper (see [40]), he proved that M k (0) < oo for all k, 
for a covariance function r having a behavior near such that 
s 2 s 3 

r(s) = 1 — — + c— + o(s 3 ) as s — ► (with c some constant > 0). 

z o . 

Later (see [23] and [44]). for X with path continuous nth derivative X^ and 
spectral distribution function F(X), he proposed a series of sufficient conditions 
involving F and o*(h) := IE[{x[^ h - X ( t n) ) 2 } for having E[{N t {0)) k ] < oo. 
Let us give an example in terms of the spectral density of X among the sufficient 
conditions he proposed. 



Theorem 2.10 (Cuzick, 1978) 

If X has a spectral density given by /(A) 

and a > 3fc/2 - 1, M t k (0) < oo. 



:=!/ (l + A 3 |logA| Q ), then for k> 2 



Those results are not restricted to zero crossings and would also apply to a large 
family of curves (see [39] and [43]). 

However, necessary conditions arc more difficult for higher monents, the main 
difficulty lying in obtaining sharp lower bounds for the of defined in (|2. 19|) . 

• Marcus (1977): generalized Rice functions. 

Marcus (see [lOlj ) generalized results of Cramer and Leadbetter and Ylvisaker 
by considering not only Gaussian processes but also by computing quantities 
such as MJSfl 1 (xi) ■ ■ ■ N'( k (x fe )] for levels , x k and integers ji,'" >Jk, 

called generalized Rice functions. 

For the proofs, the author returns to the approach used by Kac (see [73]), Ivanov 
(see [7D]) and Ito (see [7T]) to obtain the mean number of crossings at a fixed 
level, which consists in first finding a function that counts the level crossings of 
a real valued function, then in substituting X for the function and finally, in 
considering the expectation. 

• Nualart and Wschebor (1991). 

We know that the general Rice formula giving the factorial moments of the 
number of level crossings of a stochastic process satisfying some conditions can 
hold whether finite or not. 

In the search of conditions for the finiteness of moments of the number of cross- 
ings, Nualart and Wschebor (see [113] ) proposed some sufficient conditions in 
the case of a general stochastic process, that reduce in the Gaussian case to: 

Theorem 2.11 (Nualart and Wschebor, 1991) 

If X — (X t ,t G I C 1R) is a Gaussian process having C°° paths and such that 
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var(X t ) > a > 0, t G /, ifcen M t (tt) < oo /or every ieueZ it G -S? and every 
keff. 

• ^zews emd Wschebor (2001). 

The computation of the factorial moments of the crossings of a process is still 
a subject of interest. In particular, when expressing these factorial moments by 
means of Rice integral formulas (of the type of (|2.16p in the case of the 2nd 
factorial moment of the zero crossings, for instance), there arises the problem 
of describing the behavior of the integrands (appearing in these formulas) near 
the diagonal; it is still an open (and difficult) problem, even though partial 
answers have been provided, as for instance by Azai's and Wschebor (see [TO] . 
[12], and references therein). Let us give an example of the type of results they 
obtained, which helps to improve the efficiency of the numerical computation of 
the factorial moments, in spite of the restrictive conditions. 

Proposition 2.1 (Azai's and Wschebor, 2001) 

Suppose that X is a centered Gaussian process with C 2fc_1 paths (k integer > 2), 
and that for each pairwise distinct values of the parameter t\, ■ ■ • ,tk € / the 
joint distribution of (X th , X[ , • • • , AT t 2fc_1 , h = 1, 2, • • • , k) is non degenerate. 
Then, as t\,t2, ■ ■ ■ ,tk — > t* , 

/ x' 1 ---x' k Px tl ,-,x t .,x>,-,x> {0,--;0,x' 1 ---x' k )dx' 1 ---dx' k w J k {t*) JJitj-U) 4 , 
where Jk(-) is a continuous non-zero function. 



2.1.4- Two reference methods in the Gaussian extreme value theory 
• The normal comparison method. 

The main tool in the Gaussian extreme value theory, and maybe one of the 
basic important tools of the probability theory, has certainly been the so-called 
normal comparison technique. 

This method, used in the Gaussian case, bounds the difference between two 
standardized normal distribution functions by a convenient function of their 
covariances. This idea seems intuitively reasonable since the finite dimensional 
distributions of a centered stationary Gaussian process is determined by its 
covariance function. 

It was first developed by Plackett in 1954 (see |125j ). by Slepian in 1962 (see 
[Hg] ). then by Berman in 1964 and 1971 (see [22]) and by Cramer in 1967 (see 
[38] ) in the independent or midly dependent cases. An extension of this method 
to the strongly dependent case was introduced in 1975 by Mittal and Ylvisaker 
(see [109] and also the 84s paper |108j of Mittal for a review on comparison 
techniques). 

Theorem 2.12 (Normal comparison theorem; Slepian, 1962) 

Let £i = (£i(t),i G T) and £2 = {^2(t),t G T) (with T parameter set) be 

two separable normal processes (possessing continuous sample functions) having 
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same mean and same variance functions, with respective covariance function pi 
and p2- Suppose that p\{s,t) < p2(s,t), s,t€T. 



Then for any x, P 



sup£i(t) < x 



< p 



sup£ 2 (i) < x 



The proof constitutes a basis for proofs of the Bcrman inequality and a whole 
line of its generalizations, and is part of what we call the normal comparison 
method. 

Let us illustrate it by considering a pair of Gaussian vectors of dimension n, 
Xi and X 2 , that we suppose to be independent, with respective distribution 
functions Fi, i — 1,2, density functions ipi, i = 1,2, and covariance matrices 
Si = {(o'i(j,k))j j k), i = 1>2 such that cr%(j,j) — o 2 (j, ])■ Then the covariance 
matrix E/j := hY, 2 + (1 — h)Hi is positive definite. Let fh and Fh be respectively 
the n-dimensional normal density and distribution function based on S^. 
Let recall the following equation, discovered by Plackett in 1954 (see [125 ), 
recorded by Slepian in 1962 (see [146] ) and proved later in a simpler way than 
the one of the author, by Berman in 1987 (see [22]): 

Theorem 2.13 (Plackett partial differential equation, 1954) 

Let $s be the centered normal density function with covariance matrix 



{o~ij)i 



Then 



9£j 

do; 



dxidxj ' 



i +3- 



This equation will help to compute the difference between the two normal dis- 
tribution functions Fi, i = 1,2. Indeed, x being in this case a real vector with 
coordinates (xi)i<i<n, we have 



9 „ , ^ da h (i,j) 

« — 7^ F h{x) 37 ah 

da h (i,j) dh 



i>j 



fh(Vu- ■ ■■.Vi-iiXuVi+i-,- ■ ■■>yj-i,Xj,yj + t,- 



■,yn)^\dy k dh, 



r x k 

where <x> represents the integral of order n— 2: f Xl .f Xi ~ 1 f x%+1 -. .f 313-1 f Xl+ \ ..f Xn . 

T L J— oo ■/ — 00 J — 00 J — oo J — oo J — oo 

J —OO 

Thus, if for all i ^ j, a 2 (i,j) ~ o- x {i,j) > 0, then F 2 {x) - F x {x) > 0. 

The same results hold in the case of Gaussian separable functions of arbitrary 

kind. 



From those last two results, Berman obtained: 
Theorem 2.14 (Berman inequality, 1964-1992) 

Suppose that (Xu, 1 < i < n) are standard normal random variables with co- 
variance matrix Ai = (A|j) and (X 2 i,l < i < n) similarly with covariance 
matrix A 2 = (Af ■); let p^ — maxflA^ J, IAj J) and p = maxp«. Then, for any 
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real numbers x\, ■ ■ ■ ,x, 



\P [Xij ■< Xj, 1 < j < n] - P[X 2 j < x j} l < j < n]\ < (2.20) 



In particular these results hold when choosing one of the two sequences with iid 
r.v., hence the maximum does behave like that of the associated independent 
sequence; it helps to prove, under some conditions, results on the maximum 
and on the point process of exceedances of an adequate level x n of a stationary 
normal sequence with correlation function r; for instance we obtain that the 
point process of exceedances converges to a Poisson process under the weak 
dependence condition r„ logn — > (see [85], chap. 4) or to a Cox process under 
the stronger dependence condition r n logn — > 7 > 0, or even to a normal process 
if r n \ogn — > 00 (see [85], chap. 6, or [108] )). 

There is a discussion in Piterbarg (1988 for the Russian version, 1996 for the 
English one) (see |122| ) about two directions in which the Berman inequality 
can be generalized, on the one hand on arbitrary events, on the other hand for 
processes and helds in continuous time. Piterbarg points out that it is not pos- 
sible to carry the Berman inequality (|2.20|) over to the processes in continuous 
time as elegantly as it was done for the Slepian inequality (|2.12p . but provides a 
solution for Gaussian stationary processes with smooth enough paths (see the- 
orems C3 and C4, pp. 10-12 in 122]) (and also for smooth enough stationary 
Gaussian fields) . 

Finally let us mention the last refinements of the Berman inequality (|2.20|) given 

by Li and Shao in 2002 (see [89] ) that provide an upper bound in (|2.20p , cleared 
of the term (1 -p^.)-i/2. 

Theorem 2.15 (Li and Shao, 2002) 

Suppose that (Xn, 1 < i < n) are standard normal random variables with co- 
variance matrix Ai = (A,^) and (X2i,l < i < n) similarly with covariance 
matrix A 2 = (A^ ); 

let pij = max (| A i j\, | A i j | ) and p = max . Then, for any real numbers X\ , • • •, x r 



Moreover, for n > 3, for any positive real numbers x\, ■ ■ ■ , x n , and when 




< 



P[X Xj <Xj,l<j<n] 



P[X 2j < Xj, 1 < j < n] < 



(2.21) 
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assuming that > A i j > for all 1 < i, j < n, then 

P[X Xj <Xj,l<j <n] < P[X 2j < Xj,l <j < n] < (2.22) 

1 ^ /7T - 2arcsin(A 1 1 1 ) 

P[Xi j <x j ,Kj<n] ex.pl > In , / 2 J 

lij-j, f\ ^ \ Tj- — 2 arcsm(A ■) 

For other precise versions and extensions of this method, we can also refer to 
e.g. Leadbetter et al. (see [85]), Tong (see |158| ). Ledoux and Talagrand (see 
[87]) and Lifshits (see [90]). 

• The method of moments, also called the Rice method. 

This method, introduced by Rice to estimate the distribution of the maximum 
of a random signal, consists in using the first two moments of the number of 
crossings to estimate the probability of exceeding some given level by a trajec- 
tory of a (Gaussian) process, as shown in the lemma below (see [122] p. 27 and 
chap. 3). In particular it relies on the fact that the event (X < x, max X s > x) 

0<s<t 

implies the event that there is at least one upcrossing: (iV t + (x) > 1), knowing 
that the probability of more than one up/down-crossing of the level x becomes 
smaller as the level becomes larger. 

This method works only for smooth processes, but can be extended to non- 
stationary Gaussian processes (see |136] and |137| ) and to non-Gaussian pro- 
cesses. 

Let X = (X s , s £ [0,t]) be a.s. continuously differentiable with one-dimensional 
densities bounded. Then 

Lemma 2.4 

< E[N+(x)} + P[X > x] - P[ max X s > x] 

0<s<t 

< \ (E[N+(x) (N+(x) - 1)] + E[N-{x) (N-(x) - l)]) + P[X >x,X t > x], 

where N^(x) denotes the number of downcrossings of level x by X on [0,t] 
(recall that X s is said to have a downcrossing of x at sq > if for some e > 0, 
X s > x in (sq — e, s ) and X s < x in (sq, Sq + e) ). 

Indeed 

P[ max X s >x]= P[X > x] + P[X < x, max X s > x] 

0<s<t 0<s<t 

and 

P[X < x, max X s > x}=P[X Q < x, N+{x) = 1] + P[X < x, N+(x) > 2] 

0<s<i 

=P[N+{x) = l]-P[X >x,N+(x) = 1]+P[X < x,N+(x) > 2] 
=P[N+(x) = 1] + P[N+(x) > 2] - P[X a > x, N+{x) > 1], 
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with P[N+(x) = 1] = E[N+(x)} -^2klP[N+{x) = k] 

k=2 

and 1P[X > x,N+(x) > 1] < P[X a >x,X t >x] + P[N^ (x) > 2]. 

Nowwe can conclude since we also have P[N^~ (x) > 2] < IE[N^~ (x) (Nf(x) — 1)]/2 

and 

oo 

kP[N+(x) = k]- P[N+(x) > 2] < IE[N+{x) (N+{x) - l)]/2. □. 

fc=2 

Estimates previously proposed for the (factorial) moments can then be used at 
this stage of calculation. 

Recently, Azai's and Wschebor (see [13]) adapted this method to express the 
distribution of the maximum of a one-parameter stochastic process on a fixed 
interval (in particular in the Gaussian case) by means of a series (called Rice 
series) whose terms contain the factorial moments of the number of upcrossings, 
and which converges for some general classes of Gaussian processes, making the 
Rice method attractive also for numerical purpose. 

2.2. Crossings and local time 

2.2.1. Representations of the number of crossings 

Let X — (X s , s > 0) be a stationary Gaussian process defined on a probability 
space (Q,,J-,P), with mean zero, variance one and correlation function r such 
that — r"(0) = 1. We are first interested in having a representation of the number 
of x-crossings of X as a sum of multiple Wiener-Ito integrals or in terms of 
Hermite polynomials, where the nth Hermite polynomial H n can be defined as 

t 2 °° t n 
exp(tx - — ) = Y H n (x) — 

or as 

ff„(z) = (-ir e ^J_ (e - 2 /2 )5 n > . 

Let W be the standard Brownian Motion (or Wiener process). 

Let H(X) denote the space of real square integrable functionals of the process 

X. 

OO 

Recall that H(X) = 7i n , Tin being the Wiener Chaos, i.e. the closed linear 

71=0 

subspace of L 2 (ri,J r , P) generated by the random variables {H n (W(h)),h € 
L 2 (M, dx), \ \h\ \ — 1}, where W(h) is the stochastic integral of h with respect to 
W and Ho is the set of real constant functions. 

We can make use as well of the multiple Wiener-Ito integral I n defined as in 
Major (see [96]), since we have 

H n (W(h)) = I n (h® n ) VheL 2 (m,dx),s.t.\\h\\ = l. (2.23) 
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This integral operator /„ satisfies the multiplication rule, namely: 
for f p e L 2 S {W, m p ) and g q € L 2 (lR q , m"), with 

L 2 s (m n ,m n ) := {/„ e L 2 {m n ,B(M n ),m n ) : (2.24) 

/n(A) = /ri(-A), /n(A) = /„(A p (i), • • • , A p („)),Vp G SVj j , 

m™ denoting the product Borel measure on Wt n , S n the symetric group of per- 
mutations of {1, • • • , n} and A = (Ai, • • • , A„), then 



WpUM = E kl{P+ p W _ m C k p Cp p+q ^ k {f p ® k g q ), 



k=0 

where f P ®kgq denotes the average overall permutations of A-arguments of the 
function 

/p(Ai,--- , A p _/-, x±, ■ ■ ■ , Xk)g q (\ p -k+i, ■ ■ ■ , \ p + q -2k, — xi, ■ ■ ■ ,—Xk)m k (dx) 
R 

and p A q = min(p, q). 

We can introduce the Sobolev spaces JD 2 '" for a e M as in Watanabe (see |161j ). 
A functional / £ H(X) with the development 



n=l n=0 

belongs to iD 2 '" if and only if 



2 

2, a 



Y / ( l + n r\\fn\\i^\\(i-Lr /2 f\\i<^, 



oo 

where L is the operator defined on H(X) by Lf = —nl n (f n ). 

n=Q 

L coincides with the infinitesimal generator of the Ornstein-Uhlenbeck semi- 
group {T u ,u > 0} of contraction operators on H{X) defined by T u (f) — 

oo oo 

J2 e- nu I n (fn), for any / = £ /„(/„). 

n— 1 n— 

In particular we have, for [3 > 0, 

(J-L)-/»=_jf e-^~ x T u du. (2.25) 

(see [161], p.24, or [110], §1.4). 

• Slud (1991, 1994): MWI integral expansion. 

Multiple Wiener-Ito integrals (MWI) may be a tool to represent and to study 
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non-linear functionals of stationary Gaussian processes, as shown by Kalliampur 
(see [71], chap. 8). Slud first applied the stochastic calculus of MWI integral ex- 
pansions, in 1991 (see |147j ) to express the number of crossings of the mean 
level by a stationary Gaussian process within a fixed time interval [0, t], the mo- 
tivation being to obtain probabilistic limit theorems for crossing-counts, then in 
1994 (see |148j or [149 ) to extend his results to C 1 -curve crossings. 



Theorem 2.16 (Slud, 1991 and 1994) 

Let X be a mean zero, variance one, stationary Gaussian process with continuous 
spectral measure, and twice- differentiable correlation function r. 
If var {N t {ip)) < oo, then 

oo 

N t (i/>(y))= E[N t (iP(y))} + Y,I n (F n ) in L 2 (n), where 

n=l 

* in the case of a given level x (i.e. ift(y) — x, Vy), the mean of N t (x) is given 
by the Rice formula &2.1\) and 

F n (X) = e —^- J* e ls{Xl+ - +x ^ x (2.26) 
2(_r"(0))i- , ^ B _ 2I ( a; )^=^^) J2 X mi ...\ m2l ds, 

1=0 * ' Kmi<-<mai<ti 

* and in the case of a C x -curve, the mean of N t (ip) * s given by the generalized 
Rice formula |£.7| ) and 

F n (X)= / e"< Al+ " +A »>— \^?^)H n {u) - y2 l H n ^(u) (2.27) 
Jo * V 

r(-r"(0))-i_ z2y2 _ ^ 

22 A mi ---A m; / e * Hi(-zy)y l 2 dy 



l<mi<-<m[<ii 



ds, 

where A is an n-vector of coordinates Xi. 

Note that the functional N t (ip) may thus be expressed as the integral on [0, t] of 
e is(AiH hA») mu itiplied by the formal MWI expansion of the form /„(/„(., s)), 

n 

w.r.t. Lebesgue measure. 

The proof is mainly based on the discrete approximation method of Cramer 
and Leadbetter (see [38]), already used to obtain the generalized Rice for- 
mula (|2.7p . by introducing the discrete-time number of crossings A^(l, 2 _m ) 
defined in (|2.8j) . which increases to N^(l). Since by hypothesis in the theorem 
var (iV^(l)) < cxd, so is the limiting variance of iV^(l,2 -m ) (via the monotone 
convergence theorem); then, because of the orthogonal decomposition of £ 2 (f2), 
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the MWI integrands for iV^(l) are the L 2 (and a.e.) limits of the corresponding 
integrands for A^(l, 2~ m ). 

To provide the MWI integrands for 2V^(1, 2 _m ), some work on indicators of the 
type I(x >c) or I (x >a)l(x h >h), is needed, after having noted that 

I((X,-o)(Jf, +h -b)<0) = I(X s >a) + l(X s + h >b) - 21( A ' 3 >a)l(X s + ft >6)- 

The main new technical tools used for the study of indicators are a 
generalization of the Hermite polynomial expansion for the bivariate-normal 
density in (|2.28[) and the identity (|2.29|) , both enunciated below because of their 
own interest. 

Lemma 2.5 (Slud, 1994) 

Vx, y € M, k, to, n € W and |t| < 1, 



j=0 
and 



-H k+j (x)H m+j (y) = i i - e »d-' a ) (2.28) 

i +3V 7 m +^ y ; 71^2 dx k dy m 



dx dyj \l + tj \^2{l + t) 

(2.29) 

The case of a constant level is simply deduced from the general case. 
Note that Slud used a different method in 1991, when considering directly a con- 
stant level, based mainly on properties of generalized hypergeometric functions; 
indeed he expressed l((x s -x)(x s+h -x)<o) as the sum i(x s ~x<o)^(x s+h -x>o) + 
l(x s -x>o)l(x s+Jl -x<o)) obtained the MWI expansion for the indicator l(x s -x>o) 
by using first the Hermite polynomial expansion of this indicator, then by study- 
ing the asymptotical behavior of hypergeometrical functions. Then he used the 
Diagram theorem to express the products of expansions as a sum of Wiener-Ito 
integrals (see [15] or [55], p. 42; another version will also be given in terms of 
Hermite polynomials by Arcones in 1994 (see [4]), and recalled below) and fi- 
nally he used Fourier transforms for computations. 

• Kratz and Leon (1997): Hermite polynomial expansion. 

At the end of the 1990s, Kratz and Leon (see JSUJ) proposed a new and direct 
method to obtain, under some assumptions on the spectral moments of the pro- 
cess, the Hermite polynomial expansion of crossings of any level by a stationary 
Gaussian process: 

Proposition 2.2 (Kratz and Leon, 1997) 

Let X be a mean zero stationary Gaussian process, with variance one, satisfying 
- r"(0) = 1 and r (lv) (0) < 00. (2.30) 
Then the following expansion holds in L 2 (fl) 



00 I5] j-t 

^W^EEVaW^I / H q -2l{X s )H 2 l 

n — n I — PI "* 



(X s )ds, (2.31) 



9=0 1=0 
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where b k (x) := —±=e~ x2/2 H k (x). 



This approach is based on an analytical formula involving the "Dirac function" , 
which formally defines the number of crossings Nt (x) as 

N t (x) = f 5 x {X s )\X s \ds (2.32) 
Jo 

and which can be made precise when approximating the Dirac function; it makes 
then explicit formulas for MWI expansions much easier to obtain than in Slud's 
papers, mainly because expanding \X S \ in Hermite polynomials in X s rather 
than in X s quite simplify the calculations (note that at s fixed, X s and X s are 
independent). 

Note that the condition (|2.30[) . stronger than the Geman condition (|2.9p . can 
from now on be replaced by the Geman condition only because of Kratz and 
Leon's recent result (theorem 12. 7p . which makes the chosen method even more 
attractive . 

Moreover this approach is natural in the sense that formally the Dirac function 

oo 

S x has the generalized Hermite expansion 5 x (u) — bk(x)Hk(u) with bk(x) — 

k=0 

r 00 i 1 2, 1 e~ x2/2 

/ 5 x (y)-H k (y) — e-y ' 2 dy = -H k (x)^=- ; then has the 

J-oo K- \/2tT Kl V 27T 



corresponding development given by lemma [2761 below, made precise by 
approximating 5 X by <&' a x := ip^ x . 

Lemma 2.6 (Kratz and Leon, 1997) 

Let X satisfy the conditions of proposition \2. e A 

Let f S L 2 ((f>(x)dx) and let (c k , k > 0) be its Hermite coefficients. One has the 
following expansion 

.( oo oo „t 

/ f{X s )\X s \ds = Y / Y, Cka21 / Hk{X s )H 2l {X s )ds 
Jo k=o i=a Ja 

oo [g/2] f 

^2^2 c q -2ia 2 i H q -2i(X s )H2i(X s )ds, 



9=0 1=0 



where (a k , k > 0) are the Hermite coefficients of the function \.\, defined by 
2 \i/2 _/2\ 1/2 (-1)' +1 



a =(-) anda 2l =(-j > l. 



By defining ( K > L = ^2^c k a 2 i / H k (X s )H 2 i(X s )ds, we note that {C K - L )k,l 

k=0 1=0 ^° 

is a Cauchy sequence in L 2 (f2) and we deduce from the Hermite expansions of 

\x\ and f(x) that C, K ' L converges to / f(X s )\X s \ds in L 2 (Q). Then to conclude 

Jo 
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the proof of the lemma, just notice that the second expansion is a consequence 
of the orthogonality. □ 

Let us be more explicit about the proof of proposition 12. 2i by giving the main 
steps. 

We will apply the previous result (lemma |2~6|) to the function 

f:=ip aiX (y) = J2b%(x)H k (y), with b%{x) = ^— rj / e~' '^-^ H k (z)dz, 
k=o 7177 ' 

to get the Hermite expansion of N° (x) :— / (p a _ x (X s )\X s \ds, namely in L 2 (Q): 

Jo 



00 [q/2] t 

N t( x ) =EE K-2i( x ) a * / H q _ 2l (X s )H 2l (X s )ds. 

q=0 1=0 •'° 



00 [q/2] t 

Let N t (x) = b i-2i{x)a 2 i J H q - 2l (X s )H 2l (X s )ds, 

q=0 1=0 

with bk{x) := lim bZ(x) — —ip(x)Hk(x). 

<y-*o k\ 
Now we can write 

E[(N t (x) - N t (x)) 2 } < 2(E[(N t (x) - N?(x)) 2 } + E[(N?(x) - N t (x)) 2 }) 

and prove by Fatou's lemma and by Jensen inequality that E[(Nt(x)— N^(x)) 2 } — 
and via the chaos decomposition that E[{N° (x) — Nt(x)) 2 }, as a — > 0, to con- 
clude to proposition 12.21 □ 

Now by using the method of regularization of Wschebor, the result of proposition 
12.21 can be extended to a larger class of processes: 



Proposition 2.3 (Kratz and Leon, 1997). 

Let X be a mean zero stationary Gaussian process with variance one and satis- 
fying the Geman condition i2.9\) . In addition we will assume that 



\9"(t + h) - 8"{t)\ < \h\Li(h), where 6"(t) := tL(t), (2.33) 

L\(h) being an even function belonging to L ([0, S],dx). 
Then the following expansion holds 

Indeed, Wschebor's regularization method allows us to drop the strong condi- 
tion involving the fourth derivative of the covariance function of the process to 
replace it by a new smoother condition, which could be named "uniform Ge- 
man condition" , constituted by the conditions (|2.9[) and (|2.33[) together, in the 
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following way. 

By introducing the regularized process X e denned by 

■1 />oo / / \ 

XI — — J if f ) X s ds, where ip is a C 2 function with support in [—1, 1] 

we can check that proposition 12.21 applies to the number N^(x) of x-crossings 
associated to X s , then we can prove that, under the uniform Geman condition, 
N^{x) — > Nt(x) in L 2 as e — > 0, to conclude on proposition ^. 31 Note that for this 
last convergence, an important step is the use of the diagram formula that we 
already mentioned (Major, 1981) to show that the partial finite developments of 
Nf(x) 

— ^= converge to the same developments in the right hand side of (|2.34[) . 
y/-r'J(0) 

from which we deduce that for each fixed q 



E 



b q -2i(x)a 2 i 




b q -2i(x)a 2 i [ H q -2i(X s )H 2 i I , Xs , , ) ds, in probability as e — > 0. □ 
Jo \ V~ r "(0) / 



E' 

1=0 

To be more explicit, let us give a version of the Major Diagram formula in terms 
of Hermite polynomials, which provides the expectations of product of Hermite 
polynomials over a Gaussian vector, a version that can be found in Brcuer and 
Major (see [29]) or in Arcones (see [4]). 
We need to introduce some notations. 

Let G := {(j, I) : 1 < j < p,l < I < lj} be the diagram of order (Zi, • • • , l p ), 
V(G) the set of vertices (j, I) of the diagram G, r{Zi, • • • , l p } the set of diagrams 
of order (Zi, • • • , l p ), Lj = {(j, I) : 1 < I < lj} the jth level of G, {((i, I), (j, m)) : 
1 < i < j < Vi 1 < I < ^ii 1 < m < lj} the set of edges ((i, I), (j, m)) s.t. every 
vertex is on only one edge. The edges connect vertices of different levels. Given 
an edge w = ((i, I), (j, m)), let di(w) = i, d2(w) = j. 

Theorem 2.17 : Diagram formula 

Let {X\, • • • , X p ) be a Gaussian vector mean zero, variance 1 and with IE[XiXj\ = 
r(i,j), l<i,j<p- Then 



E 



E 

Ger{/i,-,z p } wgv(G) 



Yl r(di(w),d 2 (w)). 



• Let us go back to the heuristic formula (|2.32p for the number of crossings. 
Suppose that X — (X s , s > 0) is a mean zero stationary Gaussian process with 
variance one, that the function of covariance r has two derivatives and satisfies 
(12301) . 

Let = ao < ai < .... < a^-i < ctk — t be the points where the change of sign 
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of the derivative of X occurs. They are in finite number because the process has 
a finite fourth spectral moment (condition (|2.30[1 ). 
There is an x-crossing of X between on and a.i+\ if 



which implies that 



fc-i 



N t (x) = 52\l ( x ai+1>x )-*(x ai>x )l 

i=0 
fc-1 



i=0 



where Y x (u) — l[ Xi0 o)(u) is the Heaviside function whose generalized derivative 
is the "Dirac function" 5 x (u) = oo if u = x and if not. 
Therefore we can formally write 



fc-i 



Nt(x) = ]T 



i=0 



S x (X s )\X s \ds, 



which gives the formula (|2.32|) . 

Note also that (|2.32[) allows to retrieve formally some classical results, as: 
- the Rice formula: 

E[N t (x)] = / E[S X (X S )]E\X S | ds = te-* 2 ' 2 yf^O) /it ; 



Cabana's formula: 

if g is a positive function on E, if G is a primitive of g, then (see [32] ) 



E 



N t {x)g{x)dx 



E 



Sx s {x)g(x)dx 



X s I ds 



g(X s ) \X s \ds 



fc-i 



= ^|G(X Qi+1 )-G(X Qi 



i=0 



the Kac formula: 

since 5o(t) = 1, by formally applying the Fourier inversion formula, we 



have So(u) = — / cos(tu)dt, 
2tt JE 



N t (0) 



S (X S ) \X s \ds = — 

^ JEJo 



cos(£X s ) | X s | d£ds. 



M.F. Kratz/ 'Level crossings and other level Junctionals 



253 



• Kratz (2000) 

If we choose for a Gaussian process the Brownian motion B, then its number 
of ^-crossings is such that either N t (x) = oo or N t (x) — a.s., whereas its local 
time L t (x) can be defined formally by 

L t (x) = f S(x)(B s )ds, 
Jo 

and rigourously by a Gaussian approximation (see the table below) or by a uni- 

r* i 

form one, as L t (x) = lira / — li_ StS \(B s )ds a.s. (see for instance [55]). 
£ ^°Jo 2e 

More generally, let us define the local time L t (x) of a Gaussian process X as the 
density of the occupation measure of X, as for instance in Berman (see [22]). Its 
construction by limiting processes is based on the sample path properties of X , 
as first introduced by Levy in the case of the Brownian motion (see for instance 
[55]). 



We may then notice that the notions of number of crossings and local time 
present some "analogies" in their formulas (heuristic and non-heuristic) , as it is 
shown in the following table: 



Local time L x t for X = (X s , s > 0) 
(for which does exist) 


Crossings number iVf for X = (X s , s > 0) 
(for which iVf does exist) 


Sojourn time of X 
on a level x in [0,t]: 

S x (t) = [ l (jt3>x) ds:= [ Y x (X s )ds 
Jo Jo 




Local time = density of S x (i) 

Fc 

L x t = [ S x (X s )ds 
Jo 

Matl- 

the Heaviside function Y x and the 
respectively by the distribution functio 

and varianc 

Moreo\ 

[ f(X s )ds= f°° f{x)L x t dx 

Jo J -oo 


Number of crossings 

irmally 

N?= f S x (X s )\X s \ds 
Jo 

icmatically 

: "Dirac function" 5 X are approximated 

n and the density of a Gaussian r.v. mean x 

e a 2 , with a — > 0. 

er, V/ € C°°, 

pt POO 

/ f{X s )\X s \ds= / f(x)N t (x)dx 

Jo J -oo 

Banach-Kac (1925-43) formula 
(see [15] and [73]) 
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Because of this type of correspondance with the local time, we became interested 
in 'validating mathematically' the heuristic formula (|2.32j) for crossing-counts 
by looking for the appropiate Sobolev space, as done for Brownian local time 
by Nualart and Vives (see |lllj . [TT2] ). and Imkeller et al. (see [69 )). 

Those authors proved that Lf e j ^B„„' ?, < ^J* r ^ „ -, , n > where lD 2,a 

* [_ EP' a , Vp > 2, < a < 1/2 

denotes the Sobolev space over canonical Wiener space obtained via completion 

of the set of polynomials F with respect to the norms ||-F||2, a = — L) a / 2 F\\2, 

a e B. 

In what concerns the number of crossings N t (x), whereas the distribution S X (X S ) 
and the random variable | X s | belong to ID 2 ' a for any a < —1/2 and for any 
a G ]R respectively, the integral in the definition of N t (x) would appear to have 
a smoothing effect, showing N t (x) as a random variable which would belong to 
£> 2 ' Q for any a < 1/4 (see [7S] and [72]). 

2.2.2. Approximation of the local time by the number of crossings 

Several authors took an interest in the problem of approximating the local time 
of an irregular process X by the number of crossings of a regularization of this 
process. One classical regularization is the one obtained by convolution, already 
presented, defined by 

XI = rj) e * Xt, where ipe(-) = —i> (~\ and ip some smooth function. 

(2.35) 

• Wschebor (1984, 1992). 

Wschebor considered this problem for a specific Gaussian process, the Brownian 
motion (in the multiparametric case, but we will only present the one parameter 
case). 

Let W = (Wt,t > 0) be a Brownian motion and let W s be the convolution 
approximation of W defined by W £ :— W * ipei where ip e is defined in (|2.35p 
with tfj, a non-negative C°° function with compact support. 
Let N*([a,b]) be the number of crossings of level x by the regularized process 
W e on the interval [a, b]. 

Wschebor (see |162j ) showed that, for x E IR, 
Theorem 2.18 (Wschebor, 1984) 

£iljVf([a,6]) £ L(x,[a,b]), as e 0, for fc = l,2,--- (2.36) 

\W\ 2 

where L(x, [a, b]) is the Brownian motion local time at level x on [a, b]. 
Later (see [164] ). he proved the following: 



M.F. Kratz/ 'Level crossings and other level Junctionals 



255 



Theorem 2.19 (Wschebor, 1992) 

For any continuous and bounded function f , for a.e. fixed w, 
ir7TT\hr / f{x)N B (x)dx - / f(x)L{x)dx -> 0, 





as 



e 







where N £ {x) := N*{[0, 1]) and L(x) := L(x, [0, 1]). 
• Azai's, Florens-Zmirou (1987, 1990) 

Azais and Florens-Zmirou, in their 87 paper (see [8j), extended Wschebor's 
result (|2.36[1 to a large class of stationary Gaussian processes, when considering 
the L 2 convergence and the zero crossings. Under some technical conditions on 
the Gaussian stationary process X and its regularization X e := X*tp e (requiring 
a bit more than the non-derivability of X and giving some bounds on the second 
and fourth spectral moments of X E ), on its correlation function r (namely r twice 
differentiable outside a neighborhood of zero, with bounded variation at zero, 
r' and r" bounded at infinity) and on the convolution kernel if) (i.e. ip(u) and 
ip'(u) bounded by a constant time |u| -2 ), they proved that 

Theorem 2.20 (Azai's and Florens-Zmirou, 1987) 

If X admits a local time L(x, [0, T]) continuous in x at zero, then 



where \2. £ denotes the second spectral moment of X e . 

Note that Azai's in 1990 (see [5]) considered more general stochastic processes 
and provided sufficient conditions for L 2 -convergence of the number of crossings 
of some smooth approximating process X e of X (which converges in some sense 
to X) to the local time of X, after normalization. 

• Berzin, Leon, Ortega (1992, 1998), Azais and Wschebor (1996). 
Whereas Wschebor proved the a.s. convergence of the variable given in theorem 
12.191 Leon and Ortega studied the L? convergence of this variable slightly 
modified, with the rates of convergence, first in 1992 (see [88]), then in 
collaboration with Berzin, in 1998 (see [21]), under slightly different hypotheses. 
Indeed, for X = (X(t),t e [0, 1]) a centered stationary Gaussian process with 
covariance function r satisfying r{t) ~t_>o 1— C|£| 2q , < a < 1, they considered 

ib E * X 

the regularized process X e = — (where the kernel %j} e approaches the 

var(ip £ * X) 

Dirac function as e — > 0), and proved that the convergence of 



(f G L A {qb{x)dx) satisfying certain regularity conditions) , is in L 2 or is the weak 
convergence, depending on the values of a (0 < a < 1/4, 1/4 < a < 3/4 and 
3/4 < a < I). 
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Meantime, Azai's and Wschebor considered in 1996 (see [10]) the a.s. conver- 

/oo 
f(x)N e (x]dx, for any continuous function / and for X belonging 
-oo 

to a larger class of Gaussian processes than the one defined by Berzin, Leon and 
Ortega. 

Note that we did not mention all the papers on this subject (or closely related 
to it); indeed much work on the approximation of local times and occupation 
measures has been done also on specific Gaussian processes such as the fractional 
Brownian motion (we can quote for instance the work by Revesz et al. in the 
80s on invariance principles for random walks and the approximation of local 
times and occupation measures, as well as the work by Azai's and Wschebor (see 
|llj ) on the first order approximation for continuous local martingales). 



2.3. Asymptotic results 

2.3.1. The law of small numbers and rate of convergence 

• Due to the extension of the Poisson limit theorem from independent to depen- 
dent Bernoulli random variables (sec 36J), there has been the same extension in 
extreme value theory, in particular for the study of the number of (up-)crossings. 
As the level x = x(t) increases with the length t of the time interval, the up- 
crossings tend to become widely separated in time provided that there is a finite 
expected number in each interval. Under some suitable mixing property satisfied 
by the process X to ensure that the occurrences of upcrossings in widely sepa- 
rated intervals can be considered as asymptotically independent events, we get 
a limiting Poisson distribution for the number of x-upcrossings -ZV+ by using the 
Poisson theorem for dependent random variables (see [160] . [35], chap. 12, [126] . 
|116j . |117j . [2"2"] . and, for a survey and further minor improvement, [55], chap. 8 
and 9). We will give here a version of such a result, proposed in Leadbetter et 
al. (see [85], theorem 9.1.2 p. 174). 
Under the condition that 

M (x) :=^N+((0,l))]<oo, 

then N^(I) < oo a.s. for bounded /, and the upcrossings form a stationary point 
process with intensity parameter fi = p(x). The point process of upcrossings 
has properties analogous to those of exceedances in discrete parameter cases, 
namely: 

Theorem 2.21 (Leadbetter et al.) 

Let us assume that the correlation function r of the centered stationary normal 
process X satisfies 

r(s) = l- y S 2 + o(s 2 ) as s^O (2.37) 
and the Berman condition given by 



r(s)logs — > 

s— »oo 



(2.38) 
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Suppose that x and t tend to oo in a coordinated way such that 



JEN^(t) = t/i(x) — > r for some constant t > 0. 



(2.39) 



Define a time-normalized point process iV+ of upcrossings having points at s/t 
when X has an upcrossing of x at s. 

Then JV+ converges to a Poisson process with intensity t as t — > oo. 

Note that the asymptotic Poisson character of upcrossings also applies to non- 
differentiable normal processes with covariance functions satisfying 

r(r) = 1— (7|T| a +o(|r|°)j as r — > (withO < a < 2, and C positive constant), 



if we consider e-upcrossings, introduced by Pickands (see |116j ) and defined be- 
low, instead of ordinary upcrossings (for which under (|2.40p their mean number 
of any level per unit time is infinite). 

For a given e > 0, a process £ is said to have an e- upcrossing of the level x at t 
if £(s) <x, Vs e (t - e, t) and, V77 > 0, 3s e (t, t + rj), £(s) > x; 
so an e-upcrossing is always an upcrossing while the reverse is not true. 
The expectation of the number of e-upcrossings of x by £ satisfying (|2.40|) can 
be evaluated, and it can be proved that asymptotically this mean number is 
independent of the choice of e for a suitably increasing level x, which leads, 
under the Berman condition and (|2.40p . to the Poisson result of Lindgren et 
al. (1975) for the time-normalized point process of e-upcrossings (see |116j and 
[85] . chap. 12 for references and more detail). 

• Remark: Another notion related to the one of upcrossings of level x by the 
process X is the time spent over x, called also sojourn of X above x, defined in 



the table of section 2.2.1 by S x (t) = / 1,% > x \ds. An upcrossing of the level x 



marks the beginning of a sojourn above x. If there is a finite number of expected 
upcrossings in each interval, then the number of sojourns above x is the same 
as the number of upcrossings. Under some mixing conditions (recalled in the 
result of Bermangiven below) and the assumption of a finite expected number 
of upcrossings in each interval, Volkonskii & Rozanov (see 160 ), then Cramer 
& Leadbetter (see [3S]) showed, under weaker conditions, by using the reasoning 
used in the proof of the Poisson limit for the distribution of upcrossings, that 
the sojourn above X has a compound Poisson limit distribution (as the sum of 
a Poisson distributed random number of nearly independent r.v. which are the 
durations of the sojourns). As early as the 70s, Berman (see [22] for the review of 
these topics) proposed an alternative to discussing upcrossings; he introduced a 
method based on Hermite polynomial expansions to study the asymptotic form 
of the sojourns of X above a level x, on the one hand for x — > oo with fixed t, on 
the other hand for x, t — > oo in a coordinated way. He considered a larger class 
of processes (with sample functions not necessarily differentiable) , allowing a 
possible infinite expected number of upcrossings in each finite interval, to prove 
the compound Poisson limit theorem. 



(2.40) 
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More specifically, by using arguments on moments, Berman (see [22) . chap. 7) 

proved that: 

if the level function x — x(t) satisfies the asymptotic condition x(t) ~ \J 2 log t 
as t — > oo and for a covariance function r such that 1 — r is a regularly varying 
function of index a (0 < a < 2) when the time tends to 0, then 
the limiting distribution of v(t)S x {t) is a compound Poisson distribution, 
where v is an increasing positive function determined by the asymptotic form of 
1 — r{t) for t — > and the compounding distribution is uniguely defined in terms 
of the index of regular variation of 1 — r(t) for t — > 0. 

• For practical use of the asymptotic theory, it is important to know how faithful 
and accurate these Poisson approximations are, and in view of applications, we 
must study the rate of convergence carefully. It took some years before getting 
information about those rates. 

Finally, Pickersgill and Piterbarg (see [120J) established rates of order t~ v for 
point probabilities in theorem 12.211 but without giving any information on the 
size of v. 

Kratz and Rootzen (see jSJ) proposed, under the condition that r(.) and r'(.) 
decay at a specified polynomial rate (quite weak assumption, even if more re- 
strictive than r(s) logs — > 0), bounds for moments of the number of upcross- 

s — >oo 

ings and also for the rate of convergence in theorem 12.211 roughly of the order 
t~ s , with 

«J=±Ainfp( S ), (2.41) 

Z s>0 

where p(.) has been defined by Piterbarg (see |121j 1 by 

(i-K*)) 2 (0A9 , 

P[ >~ l-r(sr+r'(s)\r'(s)\- 1 ' 

Their approach proceeds on the one hand by discretization and blocking to 
return to discrete parameter cases, and on the other hand by combining the 
normal comparison method and the Stein-Chen method (see |152j and [35]) as 
developed by Barbour, Hoist and Janson (see [61], [16] and also [52]), this last 
method imposing to measure the rate of convergence with the total variation 
distance, defined between integer valued r.v.'s. X and Y by 

d{X,Y) = - ^ \H X = k )- F ( Y = k )\ = sup|F(X £ A) — P{Y e 
2 ... A 



More precisely, the authors achieved the following result on what concerns the 
rate of convergence, throughout assuming for convenience that x > 1 and t > 1: 

Theorem 2.22 (Kratz and Rootzen, 1997) 

Suppose {X(s);s > 0} is a continuous stationary normal process which satisfies 
nni\) . nuty) with 7 = 2 7 VEW\) and [21$ and that r{s) > for < s < t. 
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Then there are constants K and K' which depend on r{s) but not on x ort such 
that 



x 2+2/a 



and, for t > t > 1, 



d(N+(t),V(t»(x))) < (2.43) 



d(N+(t),V(tp(x))) < K /[ -^-L . (2.44) 



The constants K in (|2.43| and K' in (I2.44p are specified in the authors paper 
(see [84]). 

Let us present briefly the method used to prove such results. 
The discretisation consists in replacing the continuous process {X(s); s > 0} by 
a sampled version {X(jq);j = 0, 1, ...}, with q roughly chosen equal to s -1 / 2 . 
which makes extremes of the continuous time process sufficiently well approxi- 
mated by extremes of its sampled version. Nevertheless it raises a new problem: 
we may count too many exceedances. So a new parameter 9 is introduced to 
divide the interval (0,t] into blocks ((fc — 1)6, kd], k = 1, t/6 with 6 roughly 
equal to i 1 / 2 (this is what blocking method refers to). This choice of 6 makes 
the blocks long enough to ensure approximate independence of extremes over 
disjointed blocks. 

It leads to considering W the number of blocks with at least one exceedance: 
t/e 

W := V I k , where I k := J, w • w V Let A = E\W\. 

Z-^l K ' K ( II11 Kj J €((li-l)8,k8] X(jq)>x) L J 

k=l 

Then by the triangle inequality for the total variation distance d, we have 

d(N x {t),r{t^)) <d(N x {t),w) + d(w,r{x)) + d(r(\),r(t^)), (2.45) 

and we estimate the three terms in the righthand side (RHS) of (|2.45|) sepa- 
rately. 

To estimate d(N x (t), W), we introduce N x q \d) the number of exceedances of x 
by the sampled process {X(jq); jq £ (0,t]} upon an interval of length 8, and 
the essential part concerns the evaluation of the difference between N x q ^ {0) and 
the probability that we have at least one exceedance by the sampled process on 
the same interval, i.e. E[N { x q \0)\ - E[h\. To evaluate it, we use lemma [2~31 
The last term of the RHS is bounded in the same way as the first one, using 

d(V(X),V(tfi)) < |A- t(i\ (see e.g. [BT], p.12), where X — t/x= \{E\h] - On). 



To study the middle term of the RHS, we use the Stein-Chen method described 
below. It means that we have to evaluate double sums of covariances of indi- 
cators that are grouped into blocks, which implies, because of stationarity, the 
study of cov{I\, Ik),k > 2. A version of the normal comparison lemma given in 
theorem 12. 141 (by taking x = min(xi, 1 < i < n)) helps to treat the case k > 3. 
When k = 2, a new parameter 9* is introduced to ensure a quasi- independence 
between blocks. Then classical normal techniques are used. 



M.F. Kratz/ 'Level crossings and other level Junctionals 



260 



Concerning references about the approximation speed in related problems, we 
could mention the works by Borodin and Ibragimov in the 80s and 90s (see [25]). 
by Jacod (see [72]), by Perera and Wschebor (see [115) ). by Berzin and Leon 
(see [23]), etc... 

• The Stein- Chen method. 

In 1970, Stein proposed a new method to obtain central limit theorems for 
dependent r.v. (see [152) ). Chen adapted this method in 1975 to the Poisson 
approximation (see [36]). 

Let (il, B, P) be a probability space, E the expectation under P and Z : f2 — > R 
a r.v. such that 1E[Z] < oo. 

The purpose of such a method is to approximate E[Z] in order to approximate 
the c.d.f. of a real r.v. W defined on (0,0, P) by l(w<w ) '■= 
The basic approach of the Stein method consists in: 

1. determining KerE := {Y : E[Y] = 0}; 

2. searching in KerE for Z — c for some constant c; 

3. concluding E[Z] ~ c. 

Given an integer valued r.v. W, the problem is then how to determine KerE, 
i.e. to characterize the set of functions h : IN — > R such that 2E[/i(W)] = 0. To 
solve this problem, a method of "exchangeability" is applied. Let us recall the 
definition of a pair of exchangeable variables. 

Let (Jl,B,P) and (Vti,B\, P\) be two probability spaces. (X, X') is an "ex- 
changeable" pair of mappings of IP) into (fii,2?i,iPi) if 

J P[X eA] = Pi(A), \/A e B x 

\ P[x e A n X' e A'} = P[x e A' n X' e A], VA A' e B\. 

A pair (X, X') of exchangeable variables satisfies the following property: 

If T = {F : n 2 -> R, F antisymmetric (i.e. F(x, x') = -F{x' , x)); E\F(X, X')\ < oo}, 

X = {h : Q — * R measurable, s.t. E\h(X)\ < oo}, 
andT : T ^ X s.t. (TF)(X) = E[F(X,X')\X], 

then (X, X') satisfies EoT = 0, with the composition defined by (EoT)(F) := 
E[TF] . 

To apply this method to the Poisson approximation, let us recall the following 
characterization for the Poisson law: 

W is an integer valued r.v., Poisson P(X) distributed, if and only if 
for any bounded f, E[\f{W + 1) - Wf(W)] = 0. 

Hence if E[Xf{W + 1) - Wf{W)\ ~ for all bounded function / defined on IN, 
then W is nearly distributed as V(X). 

This approximative equality is often easier to check than the direct approxima- 
tion E[f(W)] ~ E[f{W)} where Z is V{\) distributed. 

Therefore, let us illustrate the Stein-Chen method when W is the number of 
occurrences of a large number of independent events. The advantage of this 
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method is that the dependent case involves only minor transformations. 

Let (Xi)i<ci< n be n independent r.v. with values in {0, 1}, with (for fixed n) 

n 

Pi := P[X, = 1], 1 < i < n, and A := ^Pi- 

i=l 

n 

Let W := Xi. We proceed in four steps: 

z=l 

1. Construction of an exchangeable pair (W, W). 

Let (X*)i<i< n be n independent r.v., independent of the Xi and s.t. 
Xi = d X*, Vl<i<n. 

Let U be a uniformly distributed r.v. in {1, • • • , n}, independent of the Xi 
and X*. 

Let W := W - X v + X£. Then (W, W) is exchangeable. 

2. Definition of an antisymmetric function F. 

Let F(X,X') := /(X')l ( x'=x +i) - f(X)l (x=x , +1) , with f : N—> M. 
Then we can apply the property of exchangeable mappings to (W, W), 
namely E [IE* [F(W, W')]] = 0, i.e. 



E[Xf(W+l)-Wf(W)] = E (f(W+l)-f(W))Y,Pi*(x j =i) 

3=1 

(2.46) 

3. Artificial solution U\h (particular case of the general method of Stein). 
Let h : IN — > Mbe bounded and let the function U\h be defined by 

U x h{w) := - K ^f-X>- W (h(j) - E[h{Z)]) , 

j=w 3 ' 

with Z Poisson V{X) distributed r.v., solution of the equation 

Xf(w + 1) - wf{w) = h(w) - E[h{Z)}. 

Taking / := U\h in (|2.46[) provides that for all bounded function h defined 
on N, 

n 

E[h(W)} = E[h(Z)] +J2p 2 j E[V x h(Wj)}, 

3=1 

where Wj := X f and V^H := [7 A ft(w + 2) - f7 A /i(w + 1)- 

j ' 

4. Estimation of d tv (£(W),V(A)). 

Let us choose h s.t. h := Iia, Ac IN, with }ia(w) := 1-i w qa)1 a bound can 
then be deduced for V\h and we obtain 



\P[W eA]-P[ZeA}\ < (lA ~) ^p; 



I) 

3=1 
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which gives the same upper bound for the distance in variation d tv (C(W) , V(X)) 
defined by d tv {C(X), C(Y)) = d tv (X, Y) := sup \P[X € A] - P[Y € A}\. 

AcM 

Note that Barbour (see [TB]) extended the Stein-Chen method by combining 
it with coupling techniques (see [140] and |141j ) to solve the general problem 
of Poisson approximations for the distribution of a sum of r.v., not necessarily 
independent, {0, l}-valued. In the case of independent indicators, Deheuvels et 
al. (see [47] and references therein) combined semigroup theory (see [140 ) with 
coupling techniques to obtain results for the Poisson approximation of sums of 
independent indicators; it has been shown to give sharper results than Barbour's 
method (see [77]). 

2.3.2. Central Limit Theorems 
• Malevich (1969), Cuzick (1976). 

In the 70s, some work was done to prove Central Limit Theorems (CLT) for the 
number of zero crossings N t (0) as t — > oo, for instance by Malevich (see [98] ) 
and by Cuzick (see [4T]). 

Cuzick gave conditions on the covariance function of X which ensure on the 
one hand a mixing condition at infinity for X and on the other hand a local 
condition for the sample paths of X; those conditions are weaker than the ones 
given by Malevich to prove the same result, although the same type of proofs is 
used. 

Theorem 2.23 (Cuzick, 1976) 
Assume that r" exists and take r(0) = 1. 
If the following conditions are satisfied: 

i) r, r" € L 2 , 

ii) the Geman condition 12. 9\) is verified, 

...... . var{N t (0)) 2 

ill) hminf — a A > 0, 

t-^oc t 

then 

(N t (0)-E[N t (0)}) d 2 

-j= — » Al (0,cr J , as t — ► oo , 

-H^ + f( aa ^r =fl -^™)*}- 

Remark: A sufficient condition to have the assumption (Hi) and which is directly 
related to the covariance structure of X can be expressed by 

The m-dependent method. 

The proof of this theorem relies on what is called the m-dependent method, 
based on an idea of Malevich, which consists in approximating the underlying 



M.F. Kratz/ 'Level crossings and other level Junctionals 



263 



Gaussian process X (and its derivatives when they exist) by an m-dependent 
process X m (i.e. such that lE[X m (s)X m (t)] = if \s — 1\ > m), in order to apply 
the CLT for the m-dependent process given in Hocffding and Robbins (see [60]). 
Note that the m-dependent process X m is obtained directly from the stochastic 
representation of X by clipping out the portion (— oo, — m/2[U]m/2, oo) of the 
integration domain and by normalizing the resulting integral. 

• Piterbarg (1978-1996). 

Applying both the method of comparison and the method of discretization, 
Piterbarg provided a central limit theorem for the number N£ (t) of upcrossings 
at level a; by a Gaussian stationary process X (see [118] or [122] ). namely 

Theorem 2.24 (CLT for the number N+(t) of upcrossings, Piterbarg, 1978) 
Let X = (X s , s € IB) be a stationary Gaussian process, mean zero, unit variance, 
with covariance function r satisfying the Geman condition 12. 9\) and 

s{\r{s)\ + \r'{s)\ + \r"{s)\)ds < oo. (2.47) 

Jo 
Then 

varN+(t) = a 2 t{\ +o(l)) as i ^ oo, 

where 



cr = yz\(p s {x,y,x,z) exp < — — — } \ dydzds 



(i 



2tt e \ 2r"(0) 



y/-r»(0)e-*°_ Q 



2vr 

Also, the central limit theorem holds for N^(t). 

For the computation of the variance of 7V+ (t) , it suffices to write 
var(N+(t)) = E[N+(t)(N+(t) -1)] + E[N+(t)] - E[(N+(t)) 2 } , and to use 
results on moments, which, combined with the condition (|2.47p . give the con- 
vergence of a 2 . 

To obtain the CLT, Piterbarg proceeded both by discretization and by 
smoothing approximation, since only a discretized in time approximating pro- 
cess does not satisfy the conditions of the CLT in discrete time, even under 

dUB . 

He then introduced a smooth enough process X S A discretized in time, in order 
to apply to it under some conditions a known result, a CLT for the number 
of upcrossings for stationary Gaussian process in a discrete time (see for in- 
stance |122j ). This approximating process X^ has been defined as X^(s) :— 

X\kA) + (s-kA) X °« k + 1 M- XS ( kA \ 

for s £ [kA, (k+ 1)A), k = 0, 1,- •• , A such that 1/A is an integer and A — > 0, 

1 f s 

and where the smoothed process X s is given by X s (s) := — / X(s + v)dv. 

o Jo 
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• Berman (1971-1992). 

Let us return briefly to the notion of sojourn time. When choosing the level 
function x = x(t) of the asymptotic order of y/2 \ogt (as t — > oo), we recalled in 
the section Law of small numbers, that the sojourn time above the level x tends 
to a compound Poisson limit, since the sojourns are relatively infrequent and 
their contributions are few but individually relatively substantial. In this case, 
the local behavior of the correlation function r was determining in the form of 
the limiting distribution. Now when choosing the level function rising at a slower 
rate (such that x(t) / \J 2 log t is bounded away from 1), the sojourns become more 
frequent and their contributions more uniform in magnitude, implying with the 

customary (in the application of CLT) normalization, namely — F : * z — 

y/var(S t (x)) 

(with IE[S t {x)] not depending on r), a normal limit distribution. Here the local 
behavior of r does play in the normalization function, but not in the form of 
the limiting distribution. 

Berman proved this CLT for two different types of mixing conditions, namely 
when the covariance function r decays sufficiently rapidly to as t — ► oo, or 
conversely at a sufficiently small rate. 

In the case of a rapid rate of decay of r, the mixing condition is based not on the 
function r itself but on the function b that appears in its spectral representation 
given by 

r(t) = b(t + s)b(s)ds, 



when supposing that the spectral distribution is absolutely continuous with 
derivative /(A), and where b is the Fourier transform in the L 2 -sense of \J /(A). 
This mixing condition is given by b £ L 1 {^L 2 , which means in particular that the 
tail of b is sufficiently small so that r tends to sufficiently rapidly. The proof 
of the CLT is based on the m-dependent method. Indeed, Berman introduced 
a family of m-dependent stationary Gaussian processes {X m (t); — oo < t < oo} 
to approach uniformly in t the original process X in the mean square sense for 
large m. Then he deduced a CLT to the normalized sojourn of X from the CLT 
of the m-dependent process (established by adapting a blocking method used 
in the proofs of CLT for dependent r.v.). 

In the case of slowly decreasing covariances, the proof of the CLT relies on a 
method specific to Gaussian processes, based on the expansion of S t {x) in a 
series of integrated Hermite polynomials and the method of moments. In fact it 
is a special case of what is known as a non-central limit theorem for Gaussian 
processes with long-range dependence (see [50], |156j and the next section Non 
central limit theorems). 

• Slud (1991, 1994) 

Introducing the theoretical tool of Multiple Wiener Ito integrals (MWI's) 
allowed some authors as Taqqu (see |154j ). Dobrushin and Major (see [50]). 
Giraitis and Surgailis (see Maruyama (see [103] ). Chambers and Slud (see 
[3"1]). ... to prove general functional central limit theorems (FCLT) (and non 



M.F. Kratz/Level crossings and other level Junctionals 



265 



central limit theorems) for MWI expansions. In the 1990s, Slud (see [147] ) ap- 
plied a general central limit theorem of Chambers and Slud (see [34 ) to provide 
Cuzick's CLT for the zero crossings N t (0) by the process X, without needing 
Cuzick's additional assumptions to get a strictly positive limiting variance. Then 
he generalized the result to constant levels (see [149] . theorem 3.1), as follows. 

Theorem 2.25 (Slud, 1994) 

Let x £ 1R be arbitrary, and let X — (X t , t > 0) be a mean 0, variance 1, station- 
ary Gaussian process with continuous spectral measure a and twice- differ entiable 

/oo 
r 2 (s)ds < oo 
-oo 

/oo 
r(s)ds < oo if x ^ 0. Then, as t — > oo, 
-oo 

(N a (f) - e-* 2 ' 2 t^?W)/n) 4 AA(0,a 2 ), 
where a 2 > is given by the expansion 

00 -1 /"CO 

a 2 = ^-I^[|/„(A)| 2 I A 1 + --- + A„ = 0] / r n (s)ds, 

n=l n ' J -°° 

with f n G L 2 (lR n ,m n ) (introduced in \2.2J$ ) defined by 



fnikn) 



[ n / 2 l II ((\Y\l-2k 

S 2fcV 2fc-l) g "- 2fc(x)g2fc( ° )tl+2fc S V ■••A„ 2fc . 

fc=0 ^ ' 1 ' l<Tii<---<n 2 fc<n 

Note that Slud extended this result to curve crossings in cases where the curve 
is periodic and the underlying process has rapidly decaying correlations (see 
[TI9"1 . theorem 6.4). 

• Kratz and Leon (2001): a general method. 

In 2001, Kratz and Leon (see [82]) introduced a general method which can be 
applied to many different cases, in particular when the dimension of the index set 
T is bigger than one, to provide a CLT as general as possible for level functionals 
of Gaussian processes X = (X t ,t £ T). This method is a combination of two 
approaches, the one developed by the authors in 1997 (see [5D]) and one derived 
from the work of Malevich (see [HE])) Cuzick (see [5T]) and mainly Berman (see 
[22] ), that consists in approaching the process X by an m-dependent process, in 
order to be able to use well-known results on m-dependent processes. Applying 
this method, a CLT is given for functionals of (X t , X t , X t , t > 0), which allows 
in particular to get immediately the CLT for the number of crossings Nt (x) of 
X, given in Slud (see [149j . theorem 3.1). 
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We suppose that the correlation function r of our (stationary Gaussian) process 
mean zero, variance one, satisfies 



r G L 1 and r {lv) G L 2 . 

Note that (12.48|) implies that r' , r" and r'" belong to L 2 as well. 

Let Z s be a r.v. independent of X s and X s (for each s fixed) such that 



= piX,. + P2Z S , with pi = 



r"(0) 

Vr<»)(0) : 



(2.48) 



,P2 = V 1 -Pi- ( 2 - 4 9) 



Note that X s and X s are independent, as well as X s and X s , which ensures the 

existence of Z s with the stated properties. 

Let be the Hermite expansion given in H(X) by 



F t X = 



oo pt / 

dqnm / H n (X s )H q -( n + m ) I 
9 =0 0<n+m< ? ^° \ 



with d qnm such that > 0, 



XI d 2 qnm n\ml(q-(n + m))\<C(q), 

0<n-\-m<q 



H m (Z s )ds, 

(2.50) 

(2.51) 



with (C{q)) q some bounded sequence. 

Besides the property of stationarity of the process and the orthogonality of the 
chaos, which will be used to simplify the computations whenever possible, let 
us give the basic points that constitute this general method. 
> Let !pQ v t(X) be the finite sum deduced from Tt{X) for q = 1 to Q, i.e. 



*j.t(*)=E E 

1 an in 



q—1 0<n+m<g 



y/t Jo 
H m (Z s ) 



I qnm (s)ds : with 



,V = » 37 (0) 

> First we show that ^(-X") can be approximated in L 2 by the r.v. Tq t (X) 
lim lim 2E[.F t (X) - T Q , t {X)] 2 = 



(2.52) 



To prove this convergence, we follow the method developed in Kratz & Leon (see 
[81] . proof of theorem 1) where two results (cited below), one of Taqqu (see [1 55] . 
lemma 3.2) and the other of Arcones (see pQ, lemma 1), helped respectively for 



the computation of expectations of the form IE 



^ dqnm,Iqnm{fy 
0<n+m<q 



and 



E 



^ ^ dqn 1 7nidqn 2 m2 Iqnimi (s)^gn 2 m2 { S ) 

0<ni+mi<q 0<ri2+m2<iJ 
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Lemma 2.7 (Arcones inequality, 1994) 

Let X = (Xi , 1 < i < d) and Y = (Yi , 1 < i < d) be two Gaussian vectors on 
]R d , mean zero, such that E[XiXj] = !E[YiYj] = Sij, 1 < i, j < d. 
Let A 1 '^ := E[XiYj) and f a function on lR d with finite second moment and 
Hermite rankr, 1 < r < oo 7 w.r.t. X. Recall that the Hermite rank of a function 
f is defined by 



rank(f) := inf < t : 3lj with ^ lj = r and E 



3 = 1 



(f(X)-E[f(X)})l[H l] (X 1 ) 



3 = 1 



Suppose * := ( sup V |r (lj) | | V ( sup Vlr^l^ < 1. Then 

E[(f(X)-E[f(X)})(f(Y)-E[f(Y)})] < ^E[(f(X)-E[f(X)}f 

In our case, yb denotes the supremun of the sum of the absolute values of the 
off-diagonal terms in the column vectors belonging to the covariance matrix for 
the Gaussian vector 

7=^(0) V^W) J 



Under the conditions on r, we have / ip (u)du < oo 



oo 

,2/ 



Lemma 2.8 ( Taqqu, 1977) 

Let p > 2 and (X\, • • • , X p ) be standard Gaussian. Then 

E[H kl (X,)-.. H kp (X p )} = \ ' ' ' n ^ ^|> = 2 ^ < fc i-A< 

otherwise 

where ^ is a sum over all indices , i q , j q € {1, 2, • • • ,p} such that 

r 

H 7^ 3li = 1, • • • ,P and there are k\ indices 1, ki indices 2, ■ ■ ■ , k p indices p. 

> Now by classical tools (the dominated convergence theorem, Fatou lemma, 
Cauchy-Schwarz inequality, ...) and Arcones inequality, we can get the limit 
variance a 2 of Tq^{X) and prove that it is finite: 

oo 

lim lim E\T Q t(X)} 2 := a 2 = V a 2 (q) < oo, where (2.53) 

Q— »oo t^oo ' ^— ' 

9=1 

r-OO 

0~ {q) ^ ^ ^ ^ dq nimi dq n2nl2 I E[Jq nimi ( i S )Lq n2rn2 ( i S)]ds. 

Q<ni+m 1 <q0<n 2 +m 2 <q 
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> A version of the m-dependent method: Berman's method. 
This method consists in approaching the process X by an m-dependent process, 
in order to be able to use well-known results on m-dependent processes. 
Define X s as a (l/e)-dependent process to approach the process X, as follows. 
Suppose X has a symmetric spectral density /; so its ltd- Wiener representation 

/oo 
e lXt (f(X)) 1/2 dW(X), where W is a complex Wiener process 
-oo 

defined on 1R. 

Let ip be defined by ip{x) = tp*(p(x), where * denotes the convolution and where 
(p is a twice differentiable even function with support contained in [—1/2, 1/2] 
(so tj) has support in [—1, 1]). We can suppose w.l.o.g. that ||v|| 2 = 1- 
By using the Fourier inversion formula we can write 

Let £ (X) = | 2 , and then introduce 



/oo 
e M (f*0 e (X))-*dW(X), 
-00 



(note that Xf for fixed t is a standard Gaussian r.v.) and its derivatives denoted 
by X [ t j)e , 0<j<2, such that 

/OO 
e iXt {iXy{f*$ e {X))?dW{X). 
-00 

So we can prove that 

E[T Q . t (X) - T Q .t(X")} 2 ->0,withe(t) -» 0, (2.54) 



■>oo 



by using Arcones inequality and some results on the correlation between the 
process (respectively its derivatives) and the (l/c")-dependent process associated 
(respectively its derivatives), obtained vhen working with the spectral represen- 
tation of the correlation functions, namely 

Proposition 2.4 (Kratz and Leon, 2001) 

(i) For all < j, k < 2, E[x[^ £ x[ k)s ] = (-l) fe 7f +fe ) (.) converge uniformly 
over compacts and in L? as e — > towards rj t k{-) '■= (— l) k r^ +k ^ (.). 

For j = k = 0, the convergence takes place in L 1 as well. 

(ii) For all < j, k < 2, rj k (.) := JE[X^ e X^] converge uniformly over com- 
pacts and in L 2 as e — > towards rj^{.). 

> Then it is enough to consider the weak convergence of the sequence TQ.t(X e ) 
towards a Gaussian r.v. as t — > 00 to get the CLT for F* . 



M.F. Kratz/ 'Level crossings and other level Junctionals 



269 



We can write 



Q 1 r t 

FlO<i.+m<, VWO 

-77 / E E d Urn^nrn(s)ds:= ^ f f Q lx^—^= 1 Z s 
V WO „ =1 n^Tl^ VWO V V- r "(°) 



g— 1 0<n+m<<? 

re M 



where is the shift operator associated to the process 
Q 

and fQ(xi,X2,X 3 ) =E E d qnm H n(xi)H q - )H m (x3). 
q—1 0<n+m<q 

Hence the weak convergence of J-q(X 6 ) towards a Gaussian r.v. is a direct 
consequence of the CLT for sums of m-dependent r.v. (see [BOj and [21]), which, 
combined with (|2~52f . (j2~53|) and (|2~54f . provides the CLT for ^(-X"), namely 

Theorem 2.26 (Kratz and Leon, 2001) 
Under the above conditions, we have 

T t (X) := ^1 ptlA _^ AA(0, a 2 ) as t -» 00, 

00 

where JE[F^} = tdaoa and a 2 = f 2 ('?) < oo? with 

q=l 

0<ni+mi<g 0<n2+m2<£ 

and (s) =H n {X 8 )H q _ 

(n+m) -, == H m (Z s ). 
\y/-r"(0)J 

Remark. Condition (12. 48|) can of course be weakened to: 

r <E L 1 and r" € L 2 (2.55) 

when considering the process X and its first derivative only, as for the number 
of crossings, 

and to r e L 1 (2.56) 

when considering the process X only. 

• As an application of theorem 12.261 we get back, under the condition (|2.55|) . 
Slud's CLT for the number of crossings of X already enunciated in theorem 
I2~25l 
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Indeed, let us consider the Hermite expansion of N t (x) given in (|2.34p . which 
corresponds to with m := 0, d qn := d gn o = b q -2l(x)a2il n =2h Vg,n, I £ IN 

and with E[F t x ] = td a0 o = tb (x)a = -e~ x2/2 . We have Vx € R, 

TT 

[9/2] 

^ b 2 q _2i{x)a\ l (2l) \(q — 2l)\ < C, with C some constant independent of q, which 

z=o 

is a consequence of proposition 3 in Imkeller et al. (see [5JJ]) given by 

Proposition 2.5 (Imkeller, Perez-Abreu and Vives, 1995) 

Let 1/4 < a < 1/2. Then there exists a constant c such that for any n G M, 

sup \H n (x)e- ax2 \ < cn~ (8Q - 1)/12 . 

Therefore we can apply theorem 12.261 to obtain 

/ N,( X) _ 

which can be written as 



\ t 7T / t^oo 



where a 2 (0 < a 2 < oo) is given by a 2 = '^^a 2 (q) 1 with 

9=1 

[9/2] [9/2] 

0-2 (^) = J! X! b q-2n 1 (x)a2n 1 b q -2n 2 (x)a2n2 X 
n!— n 2 — 



H2n 1 (Xo)H q -2n 1 ( — , ] H2n 2 {X s )Hq-2n 2 



We can easily check that a 2 ^ 0, by proving that a 2 (2) ^ (showing that the 
determinant of the matrix (expressed in the domain of frequencies) associated 
to er 2 (2) is strictly positive). We conclude to Slud's CLT. 



ds. 



2.3.3. Non central limit theorems 

From the CLT given in the previous section, or even more generally from the 
literature of the CLT for non-linear functionals of Gaussian processes, we can 
ask what happens when some of the conditions of those theorems are violated, 
in particular when we are under a condition of regular long-range dependence. 
This problem interested many authors, among whom we can cite in chronological 
order Taqqu (see |154j . [156] ). Rosenblatt (see [133] . |134j . |135j ). Dobrushin and 
Major (see [50]), Major (see [57)), Giraitis and Surgailis (see [5E]), Ho and Sun 
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(see [63]), Slud (see [148j or |149j ). MWFs expansions proved to be quite useful 
in defining the limiting behavior in non central limit theorems for functionals of 
Gaussian processes with regular long-range dependence. It is what allows Slud 
to prove, by using techniques proper to MWI's, in particular Major's non central 
limit theorem for stationary Gaussian fields with regular long-range dependence 
(see [96], theorem 8.2), the following non CLT for level-crossing counts. 

Theorem 2.27 (Slud, 1994) 

Let X be a stationary Gaussian process, mean 0, variance 1, with continuous 
spectrum and twice differentiable correlation function with regular long range 
dependence, i.e. r{s) — (1 + \s\)~ a L(s), where L is slowly varying at oo and 
< a < 1/2. Assume also that for some S € (— oo, a) and constant C < oo, for 
k > 1, 2 and all x > 0, 

d k 



\r(x)\ 



dx 



jHx) 



<ca + \x\r 



i) Let c ^ be fixed arbitrarily. Then 



N t {c)-e- c /y^(0)A d 2e- c c 2 (-r"(0)) \ 

. — > TV 0, — — r— ^ , as t — > oo. 

y/t*-°L{t) \ {1 - a){2 - a)ir 2 J ' 

ii) When c = 0. 

r 



L(t) 



«-i / . \ d -r"(0)~ / e l(Xl+x ^-l\ 



where I2 denotes the second order MWI integral operator for a stationary 
Gaussian process X with correlation function ro(s) := J e lsx o~o(dx) uniquely 
determined by 

. (1 - cosx) 2 / J 



a (dx) — / (l-x)|a; + s| a dx, s > 0. 







3. Extensions 

The aim of this section is to present examples of possible applications of some 
of the methods reviewed in the previous section. 



3.1. CLT for other non-linear functionals of stationary Gaussian 
processes 

As an application of the heuristic of the previous section, in particular of theo- 
rem [27261 we may look at various functionals related to crossing functionals of 
a stationary Gaussian process X, as for instance the sojourn time of X in some 
interval, the local time of X when it exists, or the number of maxima of the 
process in an interval. 
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• Time occupation functionals. 

The simplest application of theorem 12.261 is when the integrand appearing in 
Fj* , defined in p.50p , depends only on one variable, without needing other con- 
ditions than the smooth one r £ L 1 (see the remark following the theorem). 

> A first example of this type is when represents the sojourn of the process 
X above a level x in an interval [0, t], i.e. when := S x {t). 
It is easy to obtain the Hermite expansion of S x {t) as 



OO „f 

S x (t) =Y. d o / H i( X s)ds, 

n JO 



q=0 

1 f + °° 1 

with do — 1 — $(x) and d„ = — - / H„{u)Mu)du = -H„^i(x)4>(x), Vg > 1. 

Then an application of theorem 12. 261 yields the CLT for the sojourn time under 
the condition r £ L , namely 



r q (s)ds. 



POO x PC 

with o 2 = Y^d\\ E[H o (X )H o (X s )]ds = V q\d 2 q \ 

q=l J q=l JO 

> Another example already discussed is the local time Lf for X in the level x 

(when it exists) (see [22]). 

Its Hermite expansion is given by 

I* = ct>{x)Y,^T j H k {X s )ds = Y,h J H k {X s )ds. 



k=0 u k=0 



and again theorem 12.261 allows to retrieve its asymptotical normal behavior 
under the condition r £ L , namely 



with a 2 = Y, l l / ^[H fe (A )ff fe (A s )]d S = 2 (x) J2 77 / 
fe=i J ° fe=i *• 70 

• Number of maxima in an interval. 

One of the main concerns of extreme value theory is the study of the maximum 



r k (s)ds. 



max Xt\ of a real- valued stochastic process X — (Xt, t £ [0, TV) having con- 

fe[o,T] / 

tinuous paths, in particular the study of its distribution F. 
There is an extensive literature on this subject, going mainly in three directions, 
according to Azai's and Wschebor (see [12]): one looking for general inequalities 
for the distribution F, the other describing the behavior of F under various 
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asymptotics, and the last one studying the regularity of the distribution F. For 
more references, see [T2] . 

In the discussion of maxima of continuous processes, the upcrossings of a level 
play an important role, as was the case in the discrete case between 
maxima of sequences and exceedances of level u n through the equivalence of 
the events (M^ < u n ) = (N+ < k), M n k) being the fctb largest value of 
the r.v. X%,--- ,X n and 7V+ the point process of exceedances on (0,1] (i.e. 
N+ = #{*/n e (0, 1] : Xi > u n , 1 < i < n}). 

In the continuous case, we have already seen that crossings and maxima are 
closely related when describing the Rice method, in particular with lemma [2~H 
providing bounds of IP[ max X s > x] in terms of factorial moments of the level 

^-crossings of the process X. 

Note that Cramer (1965) noted the connection between u- upcrossings by X and 
its maximum, e.g. by {N U (T) = 0} = {M(T) < u} U {N U (T) = 0,X(0) > u}, 
which led to the determination of the asymptotic distribution of the maximum 



Recently, when looking for a reasonable way based upon natural parameters of 
the considered process X to compute the distribution of the maximum of X, 
Azai's and Wschebor (see [13]) established a method, based upon expressing the 
distribution of the maximum 2P[ max X s > x] of the process satisfying some 

regularity conditions, by means of Rice series, whose main fcth term is given 
by (— l) k+1 Vk/k\, Vk denoting the kth factorial moment of the number of up- 
crossings. This method, named 'the Rice method revisited' because inspired by 
previous works such as the one of Miroshin (see 106J ) , can be applied to a large 
class of processes, and allows a numerical computation of the distribution in 
Gaussian cases more powerful in many respects than the widely used Monte- 
Carlo method, based on the simulation of the paths of the continuous parameter 
process. 



Note also another useful connection with the maximum max X s I of a process 



Bcrman uses this equivalence between the events I max X s > x I and (S x (t) > 0) 



to study the maximum of the process X (see [22], chap. 10). 

Here we are interested in the number of local maxima of a stationary Gaussian 
process lying in some interval, and more specifically in its asymptotical behavior, 
one of the motivations being the applications to hydroscience. 



M(T). 




X, which is the sojourn time. Indeed, we can write 



max X s < x 

0<s<t 
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More precisely let ^ be the number of local maxima of X s ,0 < s < t, 
lying in the real interval [$1,^2] and let r(0) = 1. 

Kratz & Leon (1997) provided under the condition r^^(0) < 00 the Hermite 
expansion of ^ i , by adapting the proof of proposition 12.21 on the number 
of crossings, and with the change of variables (12.491) . 

Formally M£ , = I l [j3u02] {X s )8 Q (X s ) \ X s | l [0yOc) (X s )ds; more precisely: 
Jo 

Theorem 3.1 (Kratz and Leon, 1997) 
Under the condition —r^ vi '(0) < 00, we have 



<r(»)(0)^ ^ . H q _ {m+n) (0) 



1 x ' q=0 0<n+m<q v ^ v // » 



J H n (X s )H q _ {n+m) ( /_ r ^(Q) : ) Hm ( Zs ) ds 



where 5 nm is defined by 
5 nm = } 1 / / (pis + P2z)l(-oo,o)(pi^ + p 2 z)H n (x)H m (z)(j)(x)(j){z)dxdz , (3.57) 
pi and P2 satisfying \2.49\ ). 

As for the number of crossings, the condition appearing in the theorem could 
be weakened by taking a similar condition to (|2 . 33[) for the fourth derivative of 
r(.); it would then provide that IE(M^ i ^) 2 < 00. 



Previously two cases of application of theorem 12 . 261 were considered, on the one 
hand when the integrand appearing in , defined in (|2.50p . depends only on 
one variable (under the simple condition r € L 1 ), on the other hand when the 
integrand depends on two variables (under the conditions r S L 1 and r" G L 2 ), 
case of our main study, the number of crossings, but which also concerns any 
convex combination of it as, for instance, Cabana estimator (see |32j ) of the 

second spectral moment (slightly modified) given by 7 = — / N t (x)da(x) 

t J -co 

(with a(.) a distribution function on 1R), and which has been studied by Kratz 
& Leon (see [ST]). 

To apply theorem 12. 261 to obtain a CLT for the number of local maxima ^ 
made us consider the last possible case, i.e. when the integrand appearing in 
Ff depends on three variables. 

Note that the condition (0) < 00 and r € L 1 imply the condition (I2.48|) of 
theorem 12.261 and that we can write M% ^ = F* when taking 

/ rM(0) 3nmff g -(m+n)(0) . . ., , . . 

d„nm = — \ / rr ; ■ Moreover we can easily check that 

y -r"(0) [q - (m + n))V2vr 

S^Lo d qnm n ~ ( n + m )) ■ < ^ 1 w ^ n C some constant, by using again 

proposition [23] 
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Therefore theorem 12^261 yields a CLT for Mff 



L ]0i fts]' resu ^ obtained by Kratz & 
Leon in 2000 (see [81]) with another method which consisted mainly in adapting 
and verifying conditions (al)-(a3) in the 1994 Slud's paper (see |149j . p. 1362). 
But, as already stated by Berman (see [21], pp. 62-63), it wouldn't have been 
possible to benefit of this method when working in a dimension higher than one, 
which confirmed us in the idea of finding a more general method as the one put 
forward to demonstrate theorem 12.261 

Theorem 3.2 (Kratz and Leon, 2000) 

Under the conditions ~A m ^(0) < oo and r € L\, 



E[M, X 



where 



Vt 



t v/r( ra )(0) 
/2tF y/-r"(0) 



N(0,a 2 ) as t 



( ( , ( ;/)<]> 



0(A)* ( -—A 

92 



'2-k 




(3.58) 



Pi (i — 1,2) being defined in l[2.49\ ), and a 2 is given by 



„(iv) 



(0) 



V ' q=l 



X 



Pmi (Z )H n2 (X S )H, 



^g-(m 1 +n 1 )(0)-ffg-(m 2 +n 2 )(0) 

(<7 - (mi + ni))!(<7 - (m 2 + n 2 ))! 
/ X, 



j-(n2+m 2 ) I 



(is. 



This last result has applications to hydroscience, in particular in the study of 
random seas. 

Ocean waves have been studied for many years by many authors, in particular 
through a stochastic approach by M. Ochi (see [114j . for references also). In this 
approach, the randomly changing waves are considered as a stochastic process 
(so that it is possible to evaluate the statistical properties of waves through the 
frequency and probability domains) and in deep water as a Gaussian random 
process. This Gaussian property was first found by Rudnick in 1951 through 
analysis of measured data obtained in the Pacific Ocean. This process is 
assumed to be stationary. It is a common hypothesis for the study of ocean 
waves since the sea does behave in a stationary way when observed over short 
periods of time (see [33], |114j ). 

From the general CLT for the number of maxima of a stationary Gaussian 
process given in theorem l3.2i we can deduce the asymptotical behavior of func- 
tionals related to ocean waves, such as the number of local positive maxima of 
waves (studied when having a non-narrow-band-spectrum random process (see 
|114j . §3.3)), or wave amplitude, or even the amplitude associated with acceler- 
ation (see |114| . §4.3)). As an example, let us give the asymptotic behavior of 
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the number of local positive maxima of waves. 

Indeed, as shown by Ochi, consideration of the negative maxima is unnecessary 
in the maximum-point-process of waves for engineering problems, because the 
negative maxima do not contribute to the largest peak value (the extreme wave 
amplitude) expected to occur in a certain number of observations. 

M W +00) 

Let us give the asymptotic behavior of the ratio — y' , with > 0. 

M [0,+oo) 

Corollary 3.1 (Number of local positive maxima of waves, Kratz and Leon, 
2000) 

where a 2 is given by 

^ ' g=l 0<ni+mi<j 0<n 2 +m 2 <q 

(n x ( r\ n x fnw ^<?-(™i+ n i)(0)-ffg-(m 2 +«2)(0) 



+00 

E 



(q - (mi + ni))\(q - (m 2 + n 2 ))! 

X \ 1 

ds, 



y-r (0/ v/ -r "(oy 



mffiCj := E[M*(y)] 



47r v /-r"(0) 

Remark: for statistical purpose, a consistent estimator of the asymptotic vari- 
ance has also been proposed (see [SJ). 



3.2. Multidimensional case: dimension n for the index set 

In most cases, for applications in other fields, results in dimension larger than 
one are required. Generalizing results obtained in dimension one seems often 
possible in theory, but, depending upon the method, it can reveal itself more 
difficult than expected. 

There are three ways to consider the multidimensional case; in the first, the 
stochastic process is indexed by a set of dimension larger than one, in the second 
the process is such that each of its components is a random vector, and in the 
third these two previous cases are combined. 

We will be interested in the first case, on which Wschebor (sec [163J ) has been 
working and generalizing results of Cramer et al. and Marcus for a d-parameter 



M.F. Kratz/Level crossings and other level Junctionals 



277 



stochastic process X = (X t ,te M d ) having C 1 paths. 

Wschebor studied problems related to the level sets of the paths of X defined 
by '■= {t '■ X t = u}, proving at the same time some type of Rice formula. 
For more detail, see his 1985 book. 

It is in this context that the work of Kratz and Leon (see [82]) on Gaussian 
fields took place. They considered the problem of the asymptotic behavior of 
the length of a level curve of a Gaussian field, by adapting the method used in 
the one dimensional case, for random processes indexed by a set of dimension 
larger than one. 

Note that once again this study was motivated because of its various applica- 
tions, in particular for the random modelisation of the sea (see e.g. [9]). 
Consider a mean zero stationary Gaussian randon field {X s y, (s,t) G M 2 ) with 
variance one and correlation function r having partial derivatives dijr, 
for 1 < i + j < 2. 

Assume that r g L 1 satisfies r 2 (0, 0) - r 2 (s,i) ^ 0, for (s,t) ^ (0,0), and that 
do2r, d2or are both in L 2 . 

Let H(X) be the space of real square integrable functionals of the field 

(X s>t ;(s,t) e E 2 ). 

Let £q(t) ( u ) ^ e the length of {(s, t) G Q{T) : X s , t = u} the level curve at level 
u for the random field X, Q(T) being the square [-T, T] x [-T, T] and \Q(T)\ 
its Lebesgue measure. 

By theorem 3.2.5 of Federer (see [53], p. 244), we have for g G C(M) 
g(u)£% (T) (u)du = [[ g(X s , t )\\VX s , t \\dsdt. 

> JJq{t) 

We assume w.l.o.g. that A := E[d 2 X Sjt ] = JE^d^Xg^] = 1 and that X is 
isotropic, i.e. that the covariance matrix of (dioX Sit , dQiX Sit ) is of the form 

A ° N (see [1], §6.2). 



Ay 

(Note that in the non-isotropic case, i.e. when the density function of (dioX Syt , doiX S}t ) 
is A/"(0, S) with E = [ 11 12 V S^.j ^ for i ^ j, we can consider the isotropic 

process Y s t defined by Y s t = — X u v , where ( | = — XT 1 / 2 ( S 

and deduce the results for X s> t from the ones of Y s> t.) 

Under these conditions, we can obtain the chaos expansion in H(X) for £q(t) ( u ) j 
as well as a CLT for it, namely: 
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Theorem 3.3 (Kratz and Leon, 2001 ) 
In H(X), we have 



C Q(T)( U ) = X X d q-2(m,+l){u)c 2m ,2l X 

9=0 0<m+l<[q/2] 

/ / Hq-2{m+l){Xs,t)H2m{dwX s j)H 2 l(d( )1 X s j)dsdt 

J JQ(T) 

with dk(u) = -gHk(u)(p(u) (<f> being the standard normal density) and 

<w = £ £ U U) + Li/a) flEia 

pi=0p 2 =0 w ' v y v ill 

junction). 

The asymptotical behavior of £-q(t) ( u ) * s described by 

— — ^ AA(0 ' fj) ' 



wift iE[£^ (r) (u)] = IOTI ^/2) and 



a 2 = 



XX X °~q,2m 1 ,2l 1 (u)5 q: 2m 2 ,2l 2 (u) X 

9=0 0<mi+l 1 <[}/2] 0<m 2 +; 2 <[g/2] 

/ / 5?[/g,2m 1 ,2i 1 (0,0)/ (?j 2m 2 ,2i 2 (s,i)]ds(it, 

J JR 

where 5 q , 2m , 21 ( u ) d q -2( m +i){u)c 2m , 2 i, and 

Iq,2m,2l{s,t) \— H q _2(m+l) (X s ,t )H 2m {d\f)X St t )H 2 i (9()lX Sjt ) . 

The proof is an adaptation to dimension 2 of the methods used to obtain propo- 



sition 12.21 and theorem 12.261 respectively. 

Indeed we introduce Cq^(u, <j) := — J £q(t)( v ) ( I > ( dv and prove 

that it converges to £qit) ( u ) m L 2 ■ Then the generalization to dimension 2 



of lemma [2~6"1 will give us the Hermite expansion of £q(t) ( m > cr ) with coefficients 
dZ(u)c 2rn ,2i, with dZ(u) — > d k (u). 

<T— »0 

Lemma 3.1 (Kratz and Leon, 2001) 

Let f £ L 2 ((f)(u, v)dudv) and let (dj., k > 0) 6e its Hermite coefficients. 
One has the following expansion 



ff f(X s . t )\\WX s . t \\dsdt = 

J JQ(T) 
oo oo oo r r 

XXX dfeC2 ^ 2i / / H k (X s , t )H 2m (d 10 X s , t )H 2l (d Q1 X s , t )dsdt 

,,._n n ,_n J •/ OCT) 



fc=0 m=o ;=o 
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^ ^ dk-2(m+l)C2m,2l I I -Hfc-2(m+i) (Xs,t )^2m (dloX s j )H 2 l (<9oiA" Sj t )dsdt , 

fc0 0<m+K[ g /2] J Q( T ) 

where C2 m ,2i is given in theorem 



The first part of theorem 13.31 i.e. the chaos expansion follows then, exactly as 

in dimension one. 

As regards the CLT, we define 

Lq(t){u) 

and Lq, t -. as the finite sum deduced from £q(t) (u) for q = 1 to K, i.e. 

L Q(T)( U ) = mJ^ll/2 J2 S q,2m,2t( U ) [ [ I a ,2m,2l{s,t)dsdt, 

IWJI , = i < ra+ K [?/2] J JQ{T) 

We define in the same way Lq, t \(u,£) in which we consider the 1/e-dependent 
random field X s defined as 

X e att = f e «' XM \(J * ^)(A 1 ,A 2 ))W(A 1 ,A 2 ). 

We must however normalize the partial derivatives, dividing by the constant 
(JE[<9ioA^ ] 2 ) 1/2 = (i£[doiA^ ] 2 ) 1/2 = (-<9 2 or(0,0)) 1/2 to get random variables 
with variance one. 

The proof follows then the one of theorem 12.261 using our combined method. 

With respect to Rice formulas and the distribution of the maximum of Gaussian 
fields, let us also mention the recent works by Adler et al. and Taylor et al. (see 
|157j and [2]), and by Azais and Wschebor (see [14]). 



4. Conclusion 

Although this work on level crossings focuses specifically on the stationary Gaus- 
sian case, it is to be noted that much research has also been conducted on 
the non-stationary Gaussian case (notion of local stationarity, notion of non- 
stationarity but with constant variance; work on diffusions, Brownian motions, 
fractional Brownian motions, etc.), as well as on the non-Gaussian case (in 
particular when considering the class of stable processes). 

For completeness, we tried to make a large inventory of papers and books 
dealing with the subject of level crossings, even though all the references are 
not explicitcly mentioned in the synopsis above. 

Note that in case an author published a paper prior to a book, the only reference 
mentioned is that of his/her book. 
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